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CHARACTERISTICS  OF  THE  TURBULENT  BOUNDARY  LAYER 
WITH  HEAT  AND  MASS  TRANSFER  AT  M=6.7 

by 

James  E.  Danberg 

ABSTRACT:  Transpiration  from  a  porous  surface  into  a  turbulent 
boundary  layer  is  a  possible  technique  for  protecting  aero¬ 
dynamic  surfaces  from  the  severe  heating  rates  associated  with 
hypersonic  speeds.  The  effects  of  air  transpiration  and  wall 
temperature  on  skin  friction  drag,  heat  transfer,  and  other 
pertinent  boundary  layer  characteristics  were  determined  from 
measurements  on  a  porous  flat  plate  in  the  U.  S,  Naval  Ordnance 
Laboratory's  Hypersonic  Wind  Tunnel  No,  4  at  a  Mach  number  of 
6.7. 

The  effects  of  heat  and  mass  transfer  on  the  distribution  of 
Mach  number,  static  and  total  temperature,  and  velocity  were 
measured,  from  which  various  boundary  layer  parameters  were 
determined  such  as:  skin  friction  coefficient;  Stanton  number; 
total,  momentum,  disuiacement ,  and  energy  boundary  layer  thick¬ 
nesses.  Lowering  the  wall  temperature  decreased  the  total  and 
displacement  thicknesses  and  increased  the  momentum  and  energy 
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thicknesses.  All  the  boundary  layer  thicknesses  increased 
approximately  linearly  with  increasing  injection  rate. 

Velocity  distributions  in  the  fully  turbulent  layer  are 
adequately  described  by  a  Prandtl  mixing  length  analysis  in 
which  the  velocity  profile  is  transformed  into  a  form  similar  to 
the  semi-logarithmic  profile  of  incompressible  flow.  Unlike 
incompressible  flow,  the  constant  of  integration  of  the  mixing 
length  velocity  profile  is  found  to  increase  with  heat  transfer 
rate  and  decrease  with  increasing  mass  transfer  rate.  These 
effects  are  related  to  changes  in  the  laminar  sublayer  thickness. 

A  modified  Crocco  equation  adequately  describes  the  relationship 
between  the  temperature  and  velocity  distributions.  Based  on 
these  semi-empirical  equations,  the  skin  friction  was  numerically 
calculated  as  a  function  of  momentum  thickness  Reynolds  number, 
and  the  results  agreed  with  the  measurements.  The  skin  friction 
calculations  can  be  extended  to  the  evaluation  of  heat  transfer 
through  the  use  of  Reynolds'  analogy  which  was  experimentally 
verified  with  and  without  air  transpiration. 
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INTRODUCTION 

This  particular  study  is  concerned  with  the  effects  of 
mass  and  heat  transfer  on  a  hypersonic  turbulent  boundary  layer 
on  a  flat  plate.  Information  obtained  from  the  present  inves¬ 
tigation  has  possible  future  use  in  the  design  of  high  speed 
missiles,  satellite  re-entry  vehicles,  and  rocket  nozzle  cooling 
problems  because  mass  transfer  cooling  provides  a  method  for 
protecting  aerodynamic  surfaces  from  high  heating  rates. 
Transpiration  will  probably  be  most  practical  when  used  to  cool 
relatively  small  areas  of  high  heating  such  as  control  surfaces 
or  in  changing  the  thermal  radiative  or  electromagnetic  proper¬ 
ties  of  the  boundary  layer.  In  addition,  information  on  the 
transpiration  process  will  shed  light  on  the  heat  blocking 
effects  of  ablation  products.  In  any  case,  before  such  systems 
can  be  effectively  designed,  it  is  necessary  to  establish  the 
parameters  which  govern  the  phenomenon  and  the  magnitude  of 
these  parameters  associated  with  efficient  transpiration 
cooling. 

Fundamentally,  the  purpose  of  the  present  research  is  to 
provide  new  experimental  data  which  extends  the  Mach  number 
range  of  the  available  mass  transfer  data.  The  ultimate 
objective  is  to  provide  additional  evidence  from  which  the 
analysis  and  theories  can  be  checked,  improved,  or  on  which  new 
analyses  can  be  based. 
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In  approaching  this  objective,  it  is  the  intention  to 
obtain  more  than  the  overall  or  gross  effect  but  also  to  inves¬ 
tigate  how  the  effect  is  produced  in  more  detail.  In  this  way 
the  more  fundamental  assumptions  inherent  in  the  theories  can 
be  examined  and  suggestions  for  improvements  are  then  more 
apparent . 

Unfortunately,  with  the  status  of  today's  knowledge,  it  is 
not  possible  to  examine  the  effects  of  heat  and  mass  transfer 
on  the  physical  mechanism  of  turbulence.  If  the  exact  solutions 
of  the  conservation  equations  for  mass,  momentum,  and  energy 
could  be  obtained,  boundary  layer  profiles  of  velocity  and 
temperature  would  be  the  result  from  which  the  skin  friction 
and  heat  transfer  would  follow  directly.  At  the  present  time, 
however,  such  solutions  are  not  available  and  analyses  are  based 
on  generalized  velocity  and  temperature  profiles  determined  from 
experimental  data.  These  profiles  are  then  used  with  simpli¬ 
fied  and  integrated  momentum  and  energy  equations.  The  results 
are  formulas  of  heat  transfer  and  skin  friction  which  are  only 
as  good  and  general  as  the  approximate  profiles  initially 

assumed . 

The  specific  objectives  of  this  research  are: 

1.  To  obtain  new  experimental  data  on  both  skin  friction 
and  heat  transfer  with  and  without  mass  transfer  at 
a  Mach  number  of  6.7  on  a  flat  plate 
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2.  To  establish  the  Reynolds  analogy  under  the  above 
conditions 

3.  To  measure  and  analyze  the  velocity  and  temperature 
profiles  over  a  range  of  conditions  of  mass  transfer, 
heat  transfer  and  Reynolds  number 

Review  of  Theory 

Since  exact  solutions  of  the  conservation  equations  have 
not  yet  been  found,  most  analyses  are  based  on  an  approach  first 
suggested  by  L.  Prandtl.  The  analysis  is  called  the  mixing 
length  theory  which  is  in  some  ways  analogous  to  the  mean-free- 
path  concept  of  the  kinetic  theory  of  gases.  Although  the 
mixing  length  theory  has  never  been  fully  accepted,  it  has  pro¬ 
duced  quite  general  and  useful  results.  However,  before  dis¬ 
cussing  the  mixing  length  theory  some  alternate  theories  will 

be  mentioned. 

A  method  for  predicting  mass-,  momentum-,  and  energy-, 
transfer  rates  has  been  developed  by  Knuth^  in  which  a  serai- 
empirical  skin-friction  coefficient  Reynolds  number  relation 
for  a  fluid  of  constant  properties  is  derived.  The  effects  of 
mass  transfer  are  included  in  this  constant  property  solution. 
Compressibility,  heat  transfer,  and  composition  of  the  injected 
gas  are  accounted  for  by  evaluating  the  constant  property 
equations  at  a  reference  temperature  and  a  reference  composition 
derived  from  laminar  flow  results. 
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A  similar  empirical  method  has  been  suggested  by  Spalding 
in  which  a  functional  relationship  between  skin-friction  coef¬ 
ficient  and  Reynolds  number  is  assumed.  Empirical  factors  are 
introduced  which  modify  the  skin-friction  coefficient  and  the 
Reynolds  number  in  the  incompressible  friction  law  to  account 
for  compressibility,  heat  transfer,  etc.  These  factors  are 
then  considered  as  functions  of  the  boundary  layer  parameters, 
such  as  Mach  number,  wall-to-free  stream  temperature  ratio,  and 
mass  transfer.  The  functional  dependence  is  determined  from 
experimental  data. 

Although  both  of  these  methods  have  the  advantage  of  being 
simple  and  relatively  easy  to  apply,  they  are  not  very  satis¬ 
fying  from  a  theoretical  point  of  view  because  they  do  not 
relate  the  skin  friction  to  conditions  in  the  physical  boundary 
layer . 

Rott^  and  Persh'^  have  attempted  to  derive  theories  for 
mass  transfer  into  a  compressible  turbulent  boundary  layer 
employing  different  assumptions  than  those  used  in  the  mixing 
length  theory.  The  assumptions  employed  are,  however,  just  as 
arbitrary  as  those  in  the  mixing  length  theory. 

In  all  of  these  theories,  considerably  more  experimental 
data  will  be  required  before  a  positive  judgment  can  be  made 
concerning  the  various  assumptions  employed.  Because  of  the 
empiricism  in  all  theories  of  turhulence,  it  is  not  possible  to 
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extrapolate,  «ith  confidence,  beyond  conditions  where  data 
already  exists. 

Mixing  Length  Theory 

5 

Prandtl's  Inconpresslble  mixing  length  theory  and  a 

6  7 

closely  related  one  by  von  Karman  ’  have  been  extended  to 

8  9 

adiabatic  compressible  flow  by  Van  Driest  and  Wilson  .  The 

theory  was  further  extended  to  Include:  heat  transfer  by 

Harkness^^  and  Moore^^;  roughness  in  incompressible  flow  by 

12 

Prandtl  and  Schlichting  ;  roughness  In  compressible  flow  by 

13  14  15 

Fenter  ;  and  mass  transfer  by  Rubesln  ,  Dorrance  and  Dore 

16 

and  Lapin  .  Several  authors  have  extended  the  basic  theory  to 
Include  injection  of  gases  other  than  air  and  including  the 
effects  of  chemical  reactions.  For  further  discussion  on  the 
latter  theories,  see  reference  17. 

All  of  these  theories  are  based  on  a  suggestion  by  Prandtl 
that  the  turbulent  shear  stress  In  the  boundary  layer  could  be 
written  as  follows 


-  A  dii 

r  =  pv'u'  -  f  * 


ci  U 

dy 


(1) 


Essentially,  the  formulation  is  based  on  dimensional  reasoning 
and  a  very  simplified  physical  model  of  turbulence.  Prandtl 
found  from  experimental  data  that  the  mixing  length,  was 


i  =  k:/ 


(2) 
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Theoretical  arguments  proposed  by  von  Karman  lead  to  the  fol 
lowing  mixing  length  equation 


When  one  or  the  other  mixing  length  is  employed  in  equation  (1), 
and  then  equation  (1)  is  inserted  into  the  momentum  equation, 
the  result  is  a  very  complicated  differential  equation  for  the 
velocity  profile.  Since  the  equation  derived  in  this  way  can¬ 
not  be  solved,  a  major  assumption  is  made  at  this  point  in  order 
to  obtain  a  velocity  profile.  It  is  assumed  that  all  the  terms 
in  the  momentum  equation  and  the  continuity  equation  containing 
derivatives  with  respect  to  x,  the  stream  direction,  are  small 
compared  to  the  remaining  terms.  In  other  words,  the  flow  is 
reduced  to  that  of  a  Couette  flow  where  x  affects  the  local 
profile  only  through  the  wall  shear  stress.  This  assumption 
has  been  shown  to  be  untrue  except  as  an  approximation  in  a 
region  very  near  the  wall.  Nevertheless,  the  assumption  allows 
the  computation  of  a  nearly  correct  velocity  profile,  and 
although  the  x  terms  are  not  negligible  everywhere  in  the 
boundary  layer,  the  resulting  equation  is  a  reasonably  good 
approximation  except  in  the  laminar  sublayer  and  at  the  very 

outer  edge  of  the  boundary  layer. 

Compressibility  and  heat  transfer  affect  these  equations 
only  through  the  local  density  which  appears  in  equation  (1). 

It  is  on  this  point  that  the  differences  between  the  Prandtl 
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and  von  Karman  mixing  length  enter  the  problem.  In  incompress¬ 
ible  flow,  both  forms  of  the  mixing  length  give  the  same 
velocity  distribution.  When  the  density  within  the  boundary 
layer  is  not  constant,  the  resulting  velocity  distributions  are 
somewhat  different  depending  on  which  formula  is  used.  It  can 
be  shown  that  Prandtl's  equation  gives  a  slightly  higher  veloc¬ 
ity  profile  slope,  but  the  difference  is  slight.  Therefore,  it 
is  not  surprising  that  many  theories  which  are  based  on 
von  Karman 's  mixing  length  give  essentially  the  same  answer  as 
would  be  obtained  from  Prandtl’s  when  terms  of  secondary 
importance  are  dropped.  No  critical  experiment  has  been  per¬ 
formed  to  determine  which  of  the  two  formulations  is  more 
correct.  In  fact,  in  view  of  the  uncertainty  in  the  mixing 
length  theory  as  a  whole,  it  may  be  doubted  whether  a 

distinction  is  possible  at  all. 

Because  there  is  little  with  which  to  Justify  choosing 
either  equation  (2)  or  (3),  Prandtl’s  mixing  length  will  be 
used  in  analyzing  the  experimental  velocity  profiles  mainly 
because  it  is  simpler  to  work  with  when  mass  transfer  is 
involved.  The  constant  of  proportionality,  k,  will  be  taken  as 
equal  to  its  incompressible  value  of  .4.  This  assumption  is 
made  because  k  -  .4  is  well  established  in  incompressible  flow 
and  not  contradicted  by  any  compressible  data.  Exact  deter¬ 
mination  of  k  depends  on  the  derivative  of  experimental  data 
which  is  very  difficult  to  obtain  with  any  degree  of  accuracy. 
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If  the  density  and  t  are  assumed  known,  then  equation  (1) 
can  be  solved  for  the  velocity  in  the  region  of  the  boundary  layer 
where  the  theory  applies,  provided  a  constant  of  integration  is 
evaluated  from  experimental  data.  An  objective  in  analyzing 
the  present  data  will  be  to  determine  how  heat  and  mass  transfer 
affect  the  constant  of  integration. 

It  should  be  pointed  out  that  all  theories  of  transpiration 
consider  that  mass  is  transferred  uniformly  across  a  plane  which 
is  identified  as  the  surface  of  the  porous  body,  but  no  such 
porous  surface  exists.  The  porous  plate  model  developed  for  the 
present  tests  has  a  finite  number  of  pores  distributed  uniformly 
on  the  surface  so  that  mass  is  transferred  as  small  jets  of  non- 
uniform  strength  and  direction.  None  of  the  theories  have  so 
far  taken  into  account  the  destabilizing  effect  of  mass  transfer 
on  the  laminar  sublayer.  Results  from  the  present  experiments 
show  that  increased  mass  transfer  acts  on  the  sublayer  in  a  way 
similar  to  surface  roughness.  It  may  be  possible  under  some 
circumstances,  for  example,  low  Reynolds  number  and  small 
injection  rate,  that  the  skin  friction  coefficient  increases 
because  of  the  destabilizing  effect  of  blowing. 

Reynolds  Analogy 

In  the  mixing  length  theory,  it  is  necessary  to  know  the 
distribution  of  density  across  the  boundary  layer,  and  most  of 
the  analyses  which  were  cited  earlier  make  use  of  the  simple 
Crocco  integral^^  where  total  temperature  is  proportional  to 


8 


NOLTR  64-99 


velocity.  The  constants  in  the  relationship  arc  determined 
from  the  temperature  and  velocity  at  the  «all  and  outer  edKC  of 
the  boundary  layer.  The  temperature  is  related  to  density  by 
assuming  the  perfect  gas  equation  of  state  and  that  the  pressure 
is  a  constant.  Restated,  the  fundamental  assumptions  are; 

1.  At  every  point  in  the  boundary  layer  the  transport  of 
momentum  (l.e.,  the  shear  stress)  Is  proportional  to  the  flux 
of  energy. 

2.  The  constant  of  proportionality  is  the  same  at  every 

point  within  the  boundary  layer. 

3.  The  pressure  in  the  boundary  layer  is  a  constant.  The 

consequences  of  these  assumptions  are  important  because  first, 
they  imply  a  Prandtl  number  equal  to  one  everywhere,  even  in  the 
laminar  flow  at  the  wall  which  is  obviously  incorrect.  Second, 
it  is  implied  that  the  Stanton  number  is  equal  to  one  half  of 
the  skin-friction  coefficient  which  is  a  form  of  Reynolds 
analogy  for  Prandtl  number  equal  to  one. 

Insufficient  experimental  data  exists  on  which  to  base 
improvements  in  the  above  simple  theory.  Some  attempts  to 
improve  the  theories  have  been  made,  but  one  of  the  main  defi¬ 
ciencies  is  that  no  measurements  have  been  made  of  both  skin 
friction  and  heat  transfer.  Furthermore,  very  few  experiments 
in  compressible  flow  have  included  temperature  and  velocity 
profiles  and  none  exist  with  mass  transfer.  It  is,  therefore, 
the  objective  of  these  experiments  to  provide  skin  friction  and 
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heat  transier  data  at  the  same  conditions  #hich  can  be  used  to 
check  the  validity  of  Reynolds  analogy.  In  addition,  since  tem¬ 
perature  and  velocity  profiles  #ere  obtained  for  the  same 
conditions,  a  direct  investigation  of  the  Crocco  relationship 
will  also  be  made. 

Sk i n  Ft i^ ion  Law 

Once  velocity  and  temperature  profiles  have  been  determined 
theoretically,  a  skin  friction  law  can  be  calculated.  The 
integral  momentum  equation  for  the  local  skin  friction  coef¬ 
ficient  for  a  flat  plate  is 


Mixing  length  theory  applies  only  in  the  fully  turbulent  layer, 
whereas  a  laminar  region  is  generally  assumed  near  the  wall  with 
a  mixed  laminar  and  turbulent  transitional  layer  between  them. 

At  the  outer  edge,  on  the  other  hand,  there  is  another  transi¬ 
tional  region  where  the  turbulence  becomes  intermittent  and 
disappears  completely  when  the  free  stream  is  reached.  Four 
regions  are  listed  in  the  following  table  with  their  approximate 
incompressible  y  location. 
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Included  is  the  percentage  contribution  that  each  layer  makes 

3 

to  the  moroentura  thickness  at  a  Roq  =  2  x  10 

TABLE  I 

Percentage  Contribu- 

Layer  Location  tlons  t7)  Momentum 

Thickness 

Laminar  Sublayer  0  —  ^  ^  5  /U-^  0,0% 

Buffer  Layer  ^  6.8% 

Turbulent  Layer  3oS)w/u^‘^y  -  'IS  50.0% 

Velocity  Defect  Layer  ,2  S  —  V  ~  S  42.6% 

The  table  shows  that  the  turbulent  layer  and  the  defect  layer 
are  most  important,  and  it  can  be  shown  that  the  fully  turbulent 
region  increases  with  increasing  Reynolds  numbers.  The  laminar 
sublayer  and  buffer  layer  are  also  important  because  they  deter¬ 
mine  the  lower  limit  of  the  fully  turbulent  region,  but  otherwise 
they  contribute  relatively  little  to  the  momentum  thickness. 

The  velocity  defect  layer  is  physically  a  large  part  of  the 
boundary  layer  thickness  (80  percent)  but  is  approximated  rather 
well  by  extrapolating  the  fully  turbulent  profile  except  near 
the  outer  edge. 

The  principal  equations  needed  to  obtain  a  skin  friction 
law  have  been  presented  as  well  as  the  principle  assumptions 
involved.  It  can  be  concluded  that  the  fully  turbulent  layer 
is  the  single  most  important  part  of  the  boundary  layer,  and 
when  the  momentum  Integral  approach  is  used,  it  is  the  velocity 


11 


NOLTR  64-99 


and  temperature  profiles  in  this  region  which  mainly  determine 
the  skin-friction  law.  Consequently,  analysis  of  the  present 
data  is  aimed  at  determining  the  effect  of  heat  and  mass  trans¬ 
fer  on  the  fully  turbulent  layer  when  it  is  correlated  in  the 
manner  suggested  by  the  mixing  length  theory. 

Although  approximate  analytical  equations  can  be  derived 
for  the  velocity  and  temperature  profiles  within  the  assumptions 
outlined  above,  it  is  still  impossible  to  obtain  an  exact 
analytical  solution  for  the  momentum  thickness  equation  (5) 
because  of  its  complicated  integrand.  Normally  a  mathematical 
approximation  is  made  involving  a  repeated  integration  by  parts 
where  only  the  first  nonzero  term  is  retained.  The  procedure 
can  be  performed  independently  of  any  specific  profile  assump¬ 
tion  as  follows 


The  assumed  boundary  conditions  are 
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Therefore,  the  first  nonzero  term  is  equal  to 


(8) 


Van  Driest,  Rubesin,  Lapin  and  others  have  used  essentially 
the  same  type  of  approximation.  It  is  only  a  first  approxima¬ 
tion  as  can  be  seen  from  the  following  figure  which  gives  a 

geometric  interpretation  to  equation  (8)  where  the  shaded  area 

U 


is  equal  to 

^  174  ^ 


Many  of  the  details  of  the  profile  are  neglected,  but  equation 
(8)  does  show  some  of  the  basic  features  and,  therefore,  it  may 
be  expected  to  indicate  the  correct  trend  for  9  with  the  princi¬ 
pal  parameters.  It  is,  however,  only  coincidence  if  the  magni¬ 
tude  should  turn  out  to  be  correct. 

As  a  consequence  of  the  above  discussion,  it  is  apparent 
that  a  closed  form  analytical  result  employing  the  above 
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approximation  to  the  momentum  thickness  may  be  in  error  due  to 
a  number  of  factors  as,  for  example,  the  error  in  truncating 
the  expansion  or  neglecting  certain  regions  of  the  profile. 
Therefore,  no  attempt  will  be  made  to  obtain  analytical  skin- 
friction  laws  from  the  profile  correlations.  Modern  numerical 
computers  can,  of  course,  be  used  with  the  generalized  profiles 
to  obtain  practically  exact  integration  of  the  momentum  thick¬ 
ness  and  integral  momentum  equations.  Thus,  the  principle 
objectives  are  to  determine  the  effect  of  heat  and  mass  transfer 
on  the  velocity  and  temperature  profiles  and  then  to  numerically 
compute  the  skin-friction  law  based  on  generalizations  of  these 
profiles. 
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EXPERIMENTAL  FACILITY 

The  experiments  described  herein  have  been  conducted  in  the 
U.  S.  Naval  Ordnance  Laboratory’s  Hypersonic  Wind  Tunnel  No.  4, 
the  basic  components  of  which  are  described  in  reference  19, 

20,  and  21.  This  wind  tunnel  is  capable  of  operating  continu¬ 
ously  at  Mach  numbers  from  5  to  10  with  supply  temperatures  from 
300®K  to  800°K  and  supply  pressures  from  1  to  50  atmospheres. 
During  the  present  experiments,  a  Mach  number  6.7  uniform  flow, 
two-dimensional  nozzle  was  used  with  a  supply  temperature  of 
550®K  and  supply  pressures  between  15.2  and  38.0  atmospheres, 
which  is  the  maximum  possible  at  Mach  number  6.7.  Automatic 
controls  held  the  supply  conditions  constant  during  the  measure¬ 
ments  with  an  average  deviation  of  +  1®C  in  supply  temperature 
and  +0.1  atmosphere  in  supply  pressure. 

A  special  test  section  for  the  tunnel  was  used  for  the 
porous  flat  plate  experiments.  The  model  spanned  the  test 
section  from  side  wall  to  side  wall  with  its  testing  surface  on 
the  plane  of  symmetry  of  the  two-dimensional  nozzle  and  with 
the  leading  edge  located  approximately  at  the  nozzle  exit. 
Provision  was  made  in  the  test  section  wall  for  mounting  the 
traversing  mechanism  which  held  the  various  boundary  layer 
probes. 
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MODEL  DESIGN  CONSTRUCTION  AND  CALIBRATION 

The  porous  flat  plate  model  has  two  main  parts,  the  holder 
and  the  porous  insert.  The  holder  consists  of  a  sharp  leading 
edge  section  followed  by  a  flat  surface  with  provision  in  the 
center  for  the  porous  insert  and  a  plenum  chamber  underneath. 

The  porous  insert  is  designed  for  uniform  gas  injection  into 
the  boundary  layer  and  is  instrumented  for  measuring  pressures 
and  temperatures. 

Holder 

Figures  1  and  2  show  the  two  parts  of  the  holder.  Figure  2 
in  particular  shows  the  plenum  chamber  and  the  recess  for  the 
porous  insert.  On  the  right-hand  side  in  figure  1  is  the  top 
plate  which  holds  down  the  insert  after  it  is  in  place  and  pro¬ 
vides  the  sharp  leading  edge  and  continuous  smooth  (measured 
roughness  of  11+4  micro  inches)  surface  around  the  insert. 

The  leading  edge  radius  is  about  .015  mm,  and  the  underside  of 
the  plate  diverges  at  a  10®  angle  to  the  upper  surface.  Overall 
model  dimensions  are:  length,  59.06  cm;  maximum  thickness, 

3.302  cm;  and  width,  25.40  cm  at  the  leading  edge;  it  is 
wider  at  the  trailing  edge  because  the  tunnel  walls  are  tapered 
.015  cm/cm  to  account  for  the  tunnel  boundary  layer  growth. 
Figure  4  shows  a  sketch  of  the  holder  and  the  porous  insert. 

Injected  gas  is  brought  into  the  plenum  chamber  through 
the  side  wall  of  the  test  section  along  the  entire  length  of 
the  porous  insert.  Provision  is  also  made  here  for  bringing 
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out  the  thermocouple  wires  and  static  pressure  tubing  from  the 
plenum  chamber. 

Near  the  leading  and  trailing  edges,  the  holder  contains 
manifolds  for  the  porous  Insert  cooling  system.  These  manifolds 
are  soldered  to  stainless  steel  tubes  imbedded  in  the  porous 
material.  Low  temperature  oil  is  circulated  as  a  coolant 
entering  through  the  test  section  wall  directly  into  the  trail¬ 
ing  edge  manifold  and  leaving  from  the  leading  edge  manifold. 

The  back  of  the  holder  is  exposed  to  high  temperatures 
caused  by  a  complicated  shock  wave  system  behind  the  plate.  A 
cooling  system  using  circulating  water  was  Installed  on  the 
back  of  the  plate  to  prevent  the  injected  gas  in  the  plenum 
chamber  from  being  heated  and  thereby  preventing  any  conduction 
from  the  high  temperature  gas  passing  over  the  back  of  the 
holder  to  the  colder  test  surface  of  the  model. 

The  top  surface  plate  clamps  the  porous  insert  down  on 
continuous  seals  and  thereby  prevents  any  leakage  around  the 
porous  material.  As  an  added  precaution,  cement  was  put  into 
the  fine  gap  between  the  top  plate  and  the  Insert.  Measurements 
indicate  the  maximum  discontinuity  between  surfaces  before  the 
cement  was  applied  was  less  than  ,051  mm. 

Porous  Insert 

The  porous  insert  is  a  plate  48.85  cm  long,  17.78  cm  wide, 
and  1.067  cm  thick.  The  material  is  sintered  316  stainless 
steel  powder  composed  of  spheres  40  microns  and  smaller  in 
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diameter.  The  resulting  density  is  42.7  percent  of  solid  316 
stainless  steel  which  is  8.028  g/cra^.  The  testing  surface  is 
flat  within  +  .127  mm,  and  its  surface  roughness  is  60  to  70 
micro  inches. 

Temperature  of  the  insert  is  controlled  by  circulating  a 
low  viscosity  oil  through  15  equally  spaced  stainless  steel 
tubes  sintered  into  the  bottom  surface  during  manufacture  of 
the  plate.  The  axis  of  each  tube  (.343  cm  outside  diameter) 
is  parallel  with  the  long  dimension  of  the  model. 

Nine  static  pressure  orifices,  .635  mm  in  inside  diameter, 
are  located  on  the  testing  surface.  The  first  orifice  is 
located  5.08  cm  from  the  start  of  the  porous  insert  and  ,556  cm 
to  one  sif'a  of  the  plate  centerline.  The  remaining  pressure 
orifices  are  5.08  cm  apart  in  the  flow  direction  and  alternate 
on  either  side  of  the  center  coolant  tube.  Stainless  steel 
tubing  cemented  into  the  porous  material  connect  the  orifices 
with  the  read-out  system. 

Nine  thermocouple  stations  with  four  thermocouples  at  each 
station  are  located  on  the  opposite  side  of  the  center  coolant 
tube  from  the  corresponding  pressure  orifice.  The  thermocouples 
are  made  from  .254  mm  iron  and  constantan  wire.  At  a  given 
station,  the  thermocouples  are  located  as  shown  in  figure  4.  A 
view  of  the  underside  of  the  porous  insert  exposing  the  coolant 
tubes  and  manifolds  is  shown  in  figure  3. 

Thermal  conductivity — Knowledge  of  certain  properties  of  the 
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porous  material  is  necessary  in  order  to  evaluate  the  data 
obtained  from  the  porous  plate.  For  example,  one  method  to 
measure  the  heat  transfer  into  the  plate  depends  on  measuring 
the  temperature  distribution  through  the  insert  wall  and  the 
thermal  conductivity  of  the  material.  Therefore,  a  sample  of 
the  material  was  sent  to  the  National  Bureau  of  Standards  for 
thermal  conductivity  tests.  The  results  of  these  tests  are 
reported  in  reference  22,  and,  for  comparison,  the  results  of 
similar  tests  for  solid  316  stainless  steel  are  reported  in 
reference  23.  The  conductivity  of  the  porous  material  averages 
about  10  to  12.5  percent  of  solid  stainless  steel  and  has  the 
same  temperature  dependence.  At  low  temperatures  (-140®C) ,  the 
conductivity  is  2.5  x  10’^  kcal/m  sec®C  or  about  the  same  as 
glass,  and,  at  room  temperature,  the  conductivity  of  the  porous 
material  corresponds  to  that  of  granite. 

Average  mass  flow— During  the  wind  tunnel  experiments,  the  mass 
flow  through  the  porous  material  was  held  constant  at  one  of 
three  selected  values  by  regulating  the  total  flow  of  air  into 
the  system.  Although  the  mass  flow  was  held  constant,  the 
temperature  of  the  porous  material  was  varied  considerably 
during  operation,  and  this  has  a  pronounced  effect  on  the  pres¬ 
sure  drop  which  was  also  measured.  In  order  to  distinguish  the 
effects  of  changing  the  temperature  level  from  increased  contam 
ination  of  the  pores,  the  relation  between  pressure  drop,  tem¬ 
perature  distribution,  and  mass  flow  has  been  derived  and  used 
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to  reduce  all  the  data  to  common  conditions.  The  resulting 
mass  flow  per  unit  area  is: 

This  equation  is  derived  from  Darcy’s  law24.  if  the  mass  flow 
through  the  plate  is  held  constant,  but  the  temperature  is 
changed  by  varying  the  temperature  of  the  fluids  in  the  cooling 
tubes,  then  the  pressure  drop  term  decreases  as  the  average  of 
>(T  decreases.  Since  the  maximum  temperature  difference  between 
some  tests  was  as  much  as  20  percent,  the  pressure  drop  across 
the  porous  insert  changed  by  35  to  40  percent.  All  the  data 
obtained  during  calibration  and  during  wind  tunnel  tests  are 
shown  in  figure  5  where  each  value  has  been  adjusted  by  the 
above  equation  to  a  common  temperature  of  20OC.  Since  the  data 
were  accumulated  from  all  of  the  wind  tunnel  tests  and  from  the 
calibration  before  and  after  about  200  hours  of  testing,  it  is 
concluded  that  there  is  no  appreciable  uniform  contamination. 

The  calibration  tests  covered  the  same  range  of  mass  flow 
and  of  Pu^-Pw^  as  during  the  wind  tunnel  testing,  but  the  pres¬ 
sure  ratios,  p^/pu,  were  quite  different.  During  the  calibration, 
the  pressure  ratio  was  between  .936  and  1.000,  and  during  the 
wind  tunnel  tests,  it  was  between  .022  and  .033.  Despite  the 
large  difference  in  conditions,  the  data  form  a  single  curve. 

This  indicates  that  compressibility  has  very  little  effect  on 
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the  flow  through  the  pores  for  this  material. 

Local  mass  flow  calibratlon--Comparison  of  experimental  results 
with  theoretical  analyses  is  considerably  simplified  by  having 
uniform  injection  into  the  boundary  layer  and,  therefore,  the 
porous  insert  has  been  designed  to  produce  uniform  injection. 

In  order  to  verify  that  the  injection  was  actually  uniform  over 
the  surface,  the  local  mass  flow  distribution  was  measured.  In 
figure  6  the  local  measurements  on  the  plate  centerline  where 
the  boundary  layer  surveys  were  made,  and  5.08  cm  each  side  of 
the  centerline,  are  shown  as  fractions  of  the  mean  mass  flow. 
Local  mass  flow  deviates  nine  percent  from  its  mean  value.  Note 
that  the  mass  flow  is  nearly  uniform  across  the  plate. 

Calibration  tests  were  performed  before  and  again  at  the 
end  of  the  wind  tunnel  tests.  The  measurements  did  not  change 
in  form  or  magnitude.  The  maximum  deviation  in  the  second 
calibration  was  less  than  one  percent.  This  confirmed  that  very 
little  contamination  of  the  porous  material  occurred  even  after 
200  hours  of  operation. 
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INSTRUMENTATION 

The  results  presented  herein  are  based  on  boundary  layer 
surveys  made  with  Pitot  pressure  probes  and  equilibrium  tem¬ 
perature  probes.  Combined  with  surface  pressure  and  temperature 
measurements,  the  data  from  these  probes  are  sufficient  to 
describe  the  Mach  number,  total  and  static  temperature,  and 
velocity  boundary  layer  profiles. 

Pitot  Pressure  Probes 

The  Pitot  pressure  probes  used  were  made  from  circular 
stainless  steel  tubing  of  .559  mm  outside  diameter  and  an  inlet 
diameter  of  .254  mm  with  a  10°  beveled  lip.  Pressures  were 
measured  with  a  Statham  pressure  transducer  with  a  range  of 
0-1  atmosphere  (in  some  cases  a  similar  transducer  in  the  range 
0-300  mm  Hg  was  used)  and  recorded  on  a  variable  span,  strip 
chart  recorder.  The  transducer  and  recording  system  were  cali¬ 
brated  before  and  after  each  survey  against  a  mercury  manometer 
with  +  . 1  mm  measuring  accuracy. 

All  probes  used  in  the  boundary  layer  surveys  were  sup¬ 
ported  in  the  tunnel  by  a  micrometer-traversing  mechanism  which 
allowed  positioning  of  the  probe  with  an  accuracy  of  .025  mm. 
Each  traverse  was  made  from  the  free  stream  toward  the  plate, 
and  maximum  probe  movement  was  about  7.5  cm.  The  point  of  wall 
contact,  needed  to  evaluate  absolute  wall  distances,  was 
indicated  by  the  completion  of  an  electric  circuit  at  surface 
contact . 
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Several  factors  y^hich  affect  the  accuracy  of  Pitot  probe 
measurements  are  discussed  in  Appendix  I  where  it  is  concluded 
that  probe-wall  interference  is  probably  the  major  cause  of 
error  in  the  measurements  taken  near  and  at  the  point  of  wall 
contact.  Turbulent  fluctuations  also  cause  some  inaccuracies 
in  the  same  data  and,  therefore,  wall  contact  points  were 
usually  disregarded  when  describing  the  region  near  the  wall. 
Compared  to  most  of  the  data,  less  emphasis  was  given  to  the 
points  within  one  probe  diameter  of  the  wall  because  of  possible 

wall  interference. 

Pressure  Data  Reduction 

The  Rayleigh  formula^^  was  used  to  compute  the  Mach  number 
from  the  measured  Pitot  pressure  and  the  static  pressure  which 
was  assumed  constant  in  the  boundary  layer  and  equal  to  the  wall 
static  pressure.  Wall  static  pressure  was  measured  by  orifices 

in  the  plate  surface. 

The  surface  shear  stress  was  calculated  by  taking  the 
slope  of  the  Mach  number  distribution  interpolated  between  the 
probe  measurements  and  zero  Mach  number  at  the  plate  surface, 
and  multiplying  this  by  the  calculated  local  velocity  of  sound 
and  the  viscosity  of  the  air  at  the  surface  (the  viscosity  was 
obtained  from  reference  26).  Both  the  viscosity  and  the  velocity 
of  sound  just  depend  on  the  surface  temperature  of  the  plate 
which  was  indicated  by  .he  thermocouples  installed  in  the 

porous  material. 
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Temperature  Measureroent  in  Hypersonic  Boundary  Layers 

Static  temperatures  within  hypersonic  wind  tunnel  boundary 
layers  must  be  obtained  indirectly.  This  is  because  the 
presence  of  any  probe  tends  to  convert,  through  shock  waves  and 
viscosity  effects,  the  translational  energy  of  the  flow  into 
thermal  energy  in  the  immediate  vicinity  of  the  Instrument. 

Thus,  probes  tend  to  indicate  the  total  temperature  of  the  flow, 
l.e.,  the  total  energy.  Fortunately,  the  total  temperature  is 
simply  related  to  the  static  temperature  and  therefore  is  an 
adequate  measurement  when  combined  with  measurements  of  the 
flow  Mach  number. 

The  conventional  total  temperature  probe^^-Sl  jg  the 
instrument  most  frequently  adapted  for  hypersonic  boundary  layer 
work.  Such  a  probe  operates  by  adiabatlcally  compressing  a 
sample  of  the  flow  by  passing  it  through  a  normal  shock  wave. 

This  sample  enters  the  probe  and  passes  over  a  temperature 
sensing  element  (e.g.,  a  thermocouple)  at  relatively  low  speeds. 
Ideally,  the  sample  is  at  total  temperature,  and  the  sensing 
element  immersed  in  this  flow  ideally  would  indicate  total  tem¬ 
perature.  Practically,  in  a  low  density  or  low  absolute  pressure, 
high-speed  and  high-temperature  flow,  the  heat  losses  from  the 
sensing  element  due  to  conduction  and  radiation  combined  with 
the  small  heating  potential  of  the  sample  causes  large  devia¬ 
tions  in  element  temperature  from  the  ideal. 
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Equilibrium  Temperature  Probe 

In  fact,  as  the  size  of  the  total  temperature  probe  is 
decreased  or  the  flovv  density  is  decreased,  the  sensing  element 
temperature  tends  to  indicate  the  external  surface  temperature 
of  the  probe,  which  is  nearly  the  adiabatic  wall  temperature 
associated  with  the  external  boundary  layer  on  the  probe.  This 
suggests  that  by  designing  a  probe  of  suitable  geometry  where 
the  relationship  between  adiabatic  wall  temperature  and  total 
temperature  is  known,  the  measurement  of  that  adiabatic  wall 
temperature  would  be  sufficient  for  boundary  layer  work.  Such 
a  probe  might  be  called  an  adiabatic  wall  temperature  probe. 
Unfortunately,  some  conduction  and  radiation  effects  are  also 
present  in  any  practical  design  and,  therefore,  the  sensing 
element  does  not  indicate  adiabatic  wall  temperature  but  a  tem¬ 
perature  associated  with  equilibrium  between  the  heat  losses  and 
the  aerodynamic  heating.  Therefore,  such  a  probe  is  called  an 
equilibrium  temperature  probe. 

Theory  of  Equilibrium  Temperature  Probe 

A  possible  configuration  for  an  equilibrium  temperature 
probe  might  be  a  sharp,  small-angled  cone  made  from  low  emmis- 
sivity  metal  and  supported  by  a  thermal  insulator  (see  fig.  7). 
The  cone  is  mounted  in  the  flovv  with  its  axis  parallel  to  the 
flow  direction.  A  thermocouple  measures  the  temperature  of  the 
cone  (T^) •  An  additional  thermocouple  (Tb)  is  provided  in  the 
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insulating  support  to  facilitate  evaluation  of  the  conduction 
losses  effect  on  T^. 

Ideally,  the  primary  thermocouple  (T^)  indicates  the 
adiabatic  #all  temperature  (Tg)  which  is  a  property  of  the  basic 
flow,  the  cone  geometry,  and  the  type  of  boundary  layer  on  the 
cone.  Cone  adiabatic  wall  temperatures  are  directly  proportional 
to  local  total  temperature,  all  other  conditions  being  equal. 


T„  =ATe 

where 

A  =  l/fc  (I-TA,)£4(T/To)c^ 


(10) 


(11) 


(T/To)c  is  the  ratio  of  static  to  total  temperature  in  the  flow 
just  outside  the  cone  boundary  layer,  and  it  can  be  determined 
from  the  Mach  number  ahead  of  the  cone  and  the  cone  geometry 
through  use  of  conventional  cone  relations^^.  The  recovery 
factor  ( is  approximately  a  constant  because  it  depends  only 
on  the  Prandtl  number  which  varies  only  slightly  with  temperature 
under  wind  tunnel  operating  conditions.  If  the  local  cone 
Reynolds  number  does  not  exceed  the  Reynolds  number  of  transi¬ 
tion,  then  the  recovery  factor  equals  the  square  root  of  the 
Prandtl  number. 

Data  Reduction  for  Equilibrium  Temperature  Probe 

Equilibrium  temperature  probes  were  held  in  the  tunnel  by 
the  same  mechanism  as  used  for  the  Pitot  probe.  The  temperatures 
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indicated  by  the  probe  thermocouples  were  recorded  on  strip 
chan  recorders.  Wall  contact  was  indicated  by  the  conical  tip 
completing  an  electrical  circuit  when  it  first  touched  the  wall. 
Tunnel  supply  temperatures  were  recorded  during  all  surveys. 

Thus,  a  recovery  factor  for  the  probe  was  computed  for  each 
survey  when  the  probe  was  in  the  free  stream  and  assumed  constant 
for  that  survey.  Local  total  temperatures  through  the  boundary 
layer  were  calculated  by  (a)  interpolating  between  the  Mach 
number  data  to  obtain  values  corresponding  to  the  same  position 
as  the  temperature  data,  (b)  calculating  the  factor  A  in  equa¬ 
tion  (11)  from  the  cone  Mach  number  (obtained  from  ref.  25)  and 
the  constant  recovery  factor  and  (c)  calculating  Tq  from  equation 
(10). 

Heat  transfer  rates  were  calculated  from  the  slope  of  the 
total  temperature  distribution  at  the  model  surface  (equal  to 
the  slope  of  the  static  temperature  at  the  surface)  and  the 
thermal  conductivity  of  air^®  at  the  surface  temperature.  This 
procedure  supplements  the  heat  transfer  obtained  from  the  tem¬ 
pera  ure  distribution  within  the  model  wall. 

In  forming  Stanton  numbers,  the  adiabatic  wall  temperature 
is  required.  For  the  zero  mass  transfer  cases,  the  adiabatic 
wall  temperatures  were  calculated  assuming  a  recovery  factor 
equal  to  the  Prandtl  number  raised  to  the  one  third  power,  where 
wall  temperature  was  used  to  evaluate  the  Prandtl  number.  The 
effect  of  mass  transfer  on  the  recovery  factor  has  been  measured 
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by  Leadon  and  Bartle^^,  Their  value  of  recovery  factor  ratio 
(i.e.,  recovery  factor  «fith  mass  transfer  divided  by  the  zero 
mass  transfer  recovery  factor)  shown  in  figure  3  was  used  to 
obtain  the  recovery  factor  for  the  various  mass  transfer  cases. 

Static  temperature  distributions  through  the  boundary  layer 
were  calculated  from  total  temperature  and  Mach  number  distribu¬ 
tions.  Compressible  flow  tables^^  were  used  to  obtain  the 
required  temperature  ratio  at  the  corresponding  Mach  number. 

»Vith  static  temperatures  known,  local  velocities  of  sound  can  be 
calculated,  and  this  times  the  Mach  number  gives  the  local 
velocity. 
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TEST  CONDITIONS 

All  of  the  experiments  reported  herein  were  performed  in 
the  U.  S.  Naval  Ordnance  Laboratory's  Hypersonic  Tunnel  No.  4. 
Flow  Conditions 

The  Hypersonic  Tunnel  No.  4  was  equipped  with  a  fixed  block 
uniform  flow  nozzle  of  nominal  exit  Mach  number  of  6.7.  This 
nominal  exit  Mach  number  varied  slightly  with  exit  Reynolds 
number  per  unit  length.  However,  the  resulting  changes  were 
small,  as  were  the  nonuniformities  in  the  nozzle's  axial  flow. 
Both  the  effect  of  Reynolds  number  on  the  Mach  number  and  the 
uniformity  of  the  flow  are  described  in  reference  21. 

Throughout  all  tests  the  supply  temperature  was  held  con¬ 
stant  at  275®C,  while  the  supply  pressure  was  varied  between 
15  and  33  atmospheres.  The  Reynolds  number  per  unit  length 
varied  from  8  x  10®  to  19  x  10®  per  meter. 

Model  Conditions 

Four  locations  on  the  plate  were  investigated,  37.78,  42.86, 
47.94,  and  53.02  cm  from  the  leading  edge.  These  stations  were 
presumed  to  have  fully  developed  turbulent  flow  for  the  following 
reasons:  the  Mach  number  and  velocity  profiles  were  character¬ 

istic  of  supersonic  turbulent  profiles;  the  measured  wall  tem¬ 
peratures  were  characteristic  of  completed  transition,  i.e.,  a 
local  maximum  at,  or  before  station  5  (32.7  cm  from  the  leading 
idge) ;  and  profiles  obtained  nearer  the  leading  edge  than 
station  5  were  characteristically  laminar.  The  transition 
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Reynolds  number  without  air  injection  was  2.6  x  10^^.  A  rela¬ 
tively  small  decrease  was  observed  to  occur  with  air  injection, 
indicating  that  free  stream  turbulence  probably  controlled  the 
location  of  transition  more  than  any  other  factor.  Thus, 
natural  turbulent  boundary  layers  were  obtained  at  the  last  four 
stations,  and  boundary  layer  trips  were  not  required. 

The  ratio  of  wall  temperature  to  free  stream  temperature 
varied  somewhat  with  mass  transfer;  the  average  values  were 
4.1,  5.2,  and  7.6.  Because  of  the  change  in  heat  transfer  rate 
in  the  transition  region,  a  variation  in  the  wall  temperature 
occurred  in  that  region  which  could  not  be  controlled  by  the 
cooling  system.  It  was  most  pronounced  for  the  zero  mass 
transfer  and  highest  pressure  conditions  and  decreased,  along 
with  the  overall  temperature  level,  as  the  mass  transfer  rate 
increased.  Figure  9  shows  the  temperature  distribution  which 
had  the  greatest  variation.  It  amounted  to  25S:  or  about  10 
percent  of  the  driving  potential  Tg-Tw  for  heat  transfer. 

The  boundary  layer  on  the  plate  was  investigated  without 
air  injection  and  at  three  mass  transfer  rates.  These  corre¬ 
sponded  to  .09,  .17,  and  .25  percent  of  free  stream  mass  flow 
per  unit  area. 
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STATIC  PRESSURE  DISTRIBUTION 


Figure  10  shows  the  wall  static  pressures  as  measured  by 
an  oil  manometer  for  the  runs  of  15  atmosphere  supply  pressure. 
The  static  pressure  in  the  undisturbed  free  stream  was  3.65  mm  Hg 
corresponding  to  a  Mach  number  of  6.7.  It  was  observed  during 
all  the  testing  that  the  model  surface  pressure  was  considerably 
higher  than  free  stream  static  pressure.  The  amount  of  the 
excess  in  pressure  can  be  seen  from  figure  10.  In  the  down¬ 
stream  region  where  most  of  the  testing  was  performed,  the 
increase  was  due  to  the  displacement  effect  of  the  boundary 
layer.  That  is  to  say,  the  boundary  layer  was  so  thick  that  its 
effect  on  the  external  flow  could  not  be  neglected. 

In  the  zero  injection  case,  there  was  a  slight  favorable 
pressure  gradient  in  the  flow  direction  of  -4  x  lO"*^  mm  Hg/mm. 

The  effect  produced  on  the  skin  friction  coefficients  compared 
with  the  ideal  flat  plate  can  be  calculated  from  the  integral 
momentum  equation  which  can  be  written^^  as 


0  9 


(12) 


The  last  term  in  equation  (12)  represents  the  pressure  gradient 
effect.  The  second  term  on  the  right  is  approximately  the 
difference  between  the  friction  relation  for  the  flat  plate  and 
the  two-dimensional  body  which  is 

__  d 

^  dX  (13) 
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A  value  computed  for  a  typical  station  is: 

P  U|i  d  X  ^  ®  ' 

Vs 

or  about  1.5  percent  of  the  average  skin  friction  coefficient 
(cf  =  14  X  10"4).  Although  the  wall  shear  stress  obtained  from 
the  Mach  number  curves  should  be  corrected  for  the  pressure 
gradient  effect,  the  inaccuracy  in  the  determination  of  Cf  is 
much  larger  than  1.5  percent  (10  to  15  percent  based  on  repeat¬ 
ability  of  the  measurements)  and,  hence,  no  correction  was  made 

in  the  data. 

Free  Stream  Static  Pressure 

In  figure  10  for  the  mass  transfer  cases,  wall  static 

pressures  are  seen  to  rise  sharply  near  the  start  of  the  porous 
material.  The  high  pressure  in  this  region  may  be  attributed 
to  a  rapid  change  in  displacement  thickness  where  the  injected 
air  first  influences  the  boundary  layer.  Changes  in  displace¬ 
ment  thickness  are  so  rapid  that  a  shock  wave  is  formed  and 
propagates  downstream  essentially  parallel  to  the  leading  edge 
shock  wave.  It  can  be  clearly  seen  in  figure  12  where  static 
pressure  probe  surveys  perpendicular  to  the  plate  are  presented. 
The  static  pressure  surveys  were  made  for  each  station  and  with 
and  without  mass  transfer  using  a  1.14  mm  diameter  static  pres¬ 
sure  probe  with  the  orifices  7.62  cm  behind  the  sharp  tip. 

In  the  zero  injection  case,  there  is  a  distinct  peak  in  each 
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survey  *hlch  can  be  correlated  .ith  the  dovnstream  propagation 
ol  a  weak  leading  edge  shock  wave.  In  the  mass  transfer  cases 
there  are  two  peaks.  The  first  Is  the  leading  edge  shock  wave, 
and  the  second  Is  a  stronger  wave  emanating  from  near  the  start 
of  the  porous  section.  It  should  be  noted  that  whenever  the 
probe  penetrates  the  shock  wave,  there  will  be  shock-induced 
boundary  layer  separation  on  the  probe,  and  the  pressure  readings 
are  not  correct  local  static  pressures.  Nevertheless,  the 
maximum  pressures  are  associated  with  the  passage  ol  the  shock 
wave  over  the  probe  orifices.  Thus,  the  maximum  pressure  can 

be  used  to  locate  the  shock  waves. 

The  fact  that  the  origin  of  the  mass  transfer  induced 

shock  wave  is  displaced  from  the  beginning  of  the  porous 
material  is  probably  due  to  the  flexibility  of  the  long  and 
slender  static  pressure  probe.  Unsyraraetr ical  loading  on  the 
probe  tended  to  bend  it  about  5  mm  away  from  the  wall,  which 
accounts  for  the  probe  location  data  being  apparently  5  ram  too 
near  the  wall.  The  same  effect  shows  up  in  locating  the  leading 
edge  shock  wave  which  appears  more  oblique  than  it  should  a 

free  stream  Mach  number  of  6.7, 

The  flat  plate  surface  Is  shown  full  scale  In  figure  12  and 

at  each  station  on  the  plate  Is  shown  the  static  pressure 
measurement  divided  by  the  free  stream  static  pressure.  This 
pressure  ratio  Is  plotted  against  distance  from  the  wall 
which  is  also  shown  full  scale.  The  total  boundary  thickness  Is 
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indicated  by  a  dashed  line  for  the  last  four  stations  where 
profile  surveys  were  made.  An  expansion  wave  traceable  to 
imperfect  continuation  of  the  nozzle  contour  into  the  test 
section  can  be  seen  interacting  with  the  flow  at  the  most  rear¬ 
ward  stations.  The  expansion  wave,  as  can  be  seen  in  the  figure, 
did  not  reach  the  boundary  layer  in  these  tests. 

It  is  interesting  to  note  that  the  measurements  show  a 
small  pressure  rise  near  the  outer  edge  of  the  boundary  layer. 

The  higher  pressure  corresponds  quite  well  with  the  wall  static 
pressure  and,  therefore,  supports  the  assumption  that  the  static 
pressure  is  constant  across  the  boundary  layer.  These  results 
should  be  considered  qualitative,  however,  because  the  forces  on 
the  flexible  probe  may  have  caused  some  unwanted  deflections. 

In  addition,  the  point  where  the  shock  waves  reflected  from  the 
wall  begin  to  interfere  with  the  measurements  was  not  determined 
although  very  near  the  wall  (on  the  order  of  the  probe  diameter), 
it  was  possible  to  see  some  interference  effects  in  the  data. 
Nevertheless,  there  is  a  measurable  increase  in  pressure  near 
the  boundary  layer  edge  which  might  be  explained  by  the  effect 
of  turbulence  level  on  static  pressure  measurements.  The  pres¬ 
sure  rise  might  also  be  related  to  the  turning  of  the  mean 
streamlines  as  they  enter  and  penetrate  into  the  boundary  layer. 
At  the  high  Mach  numbers  of  the  outer  layers,  the  turning  could 
account  for  the  observed  pressure  increase.  In  any  case,  for 
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the  purpose  of  data  reduction,  the  static  pressure  was  considered 
constant  across  the  boundary  layer  and  equal  to  the  wall  pres¬ 
sure. 

Following  the  high  pressure  region  near  the  start  of  the 
porous  material,  there  is  a  minimum  in  pressure  which  can  be 
seen  in  figure  10  in  the  mass  transfer  cases.  The  small  pres¬ 
sure  rise  after  the  minimum  is  associated  with  the  transition 
of  the  boundary  layer  from  laminar  to  turbulent  flow.  A  similar 
effect  has  been  observed  by  Coles^S  in  connection  with  supersonic 
turbulent  boundary  layer  measurements.  The  pressure  change  is 
slight  for  the  15  atmosphere  zero  mass  transfer  cases  shown  in 
figure  10.  Figure  11  shows  a  38  atmosphere  zero  mass  transfer 
case  on  an  expanded  scale.  The  transition  region  now  extends 
from  station  3  to  6  as  indicated  by  wall  temperature  measure¬ 
ments.  With  mass  transfer  the  boundary  layer  thickness  is  con¬ 
siderably  larger,  and  the  effect  of  transition  produces  a 
measurable  pressure  change. 

Analysis  of  Displacement  Thickness  Effect 

The  displacement  effect  on  surface  pressures  in  the  region 
where  measurements  were  made  can  be  analyzed  in  the  following 
way.  Consider  a  control  volume  of  unit  width,  one  side  being 
the  plate  surface  extending  from  the  leading  edge  to  a  given 
station,  a  second  side  is  identical  and  parallel  to  the  first 
side  but  displaced  away  from  the  surface  sufficiently  to  be 
entirely  outside  the  boundary  layer.  All  the  mass  that  enters 
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the  upstream  face  of  the  control  volume  must  be  accounted  for 
as  it  leaves  the  volume.  The  total  inflow  cannot  leave  the 
downstream  face  of  the  volume  since  the  effect  of  the  boundary 
layer  is  to  decrease  the  mass  flow  near  the  plate.  Thus,  the 
flow  is  displaced  upward  away  from  the  plate;  that  is,  it  is 
given  a  component  of  velocity  away  from  the  wall  and  leaves  the 
control  volume  through  the  upper  surface.  At  supersonic  speeds, 
a  turning  of  the  flow  is  accompanied  by  a  change  in  pressure. 

As  the  flow  turns  in  the  boundary  layer,  it  is  compressed  as  in 
the  case  of  flow  through  an  oblique  shock  wave,  and  the  larger 
the  turning  angle  the  higher  the  downstream  pressure  becomes. 

At  hypersonic  speeds  only  a  very  small  turning  angle  is  required 
to  produce  measurable  amounts  of  compression. 

Qualitatively,  the  displacement  effect  can  be  calculated 
from  the  following  analyses.  The  continuity  equation  is: 


(15) 


This  equation  is  integrated  from  the  wall  to  the  edge  of  the 
boundary  layer,  y  =  6,  to  obtain 


PV 


y-s 

y=o 


(16) 


fhe 

y  = 


limits 

0.  It 


can  be  written 
can  be  assumed  that 

s  s 


at  y  =•  (5  and 
is  constant. 


pv= 

so  that : 
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fee.".-*'' 


(17) 


Therefore,  equation  (16)  can  be  written  as: 

s 


(18) 


If  equation  (18)  is  divided  by  p^UjOne  obtains 


dS' 

7  ~  ^  dx 


ft  _ 


(19) 


where 


=  Pvy^  / P« 


s*=  l'-a.Vy 


(20) 


(21) 


0 

if  the  angle  of  flow  deflection  due  to  the  displacement  effect 
of  the  boundary  layer  is  equal  to  a,  then 


v. 


u. 


=  ^'SY\  flC 


(22) 


Equation  (22)  and  (19)  show  that  the  flow  at  the  boundary  layer 
edge  will  turn  in  direct  proportion  to  the  mass  transfer  coef¬ 
ficient  and  the  change  in  slope  of  the  displacement  thickness. 

The  deflection  of  the  flow  is  associated  with  a  pressure  rise, 
and  the  applicable  relation  is  the  same  as  that  used  to  calculate 


deflection  angles  of  oblique  shock  waves. 
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This  relation  is  more  clearly  understood  if  equation  (23)  is 
approximated  by  the  more  simple  relation; 

tatiK  (24) 

which  holds  for  small  pressure  changes  00  Com- 

o 

bining  equations  (19)  and  (24),  the  final  result  is: 

y‘=  (25) 

'•0 

Equation  (25)  shows  that  the  pressure  rise  is  of  the  order  of 
t*M^times  the  turning  angle.  For  the  present  investigation, 
the  percent  pressure  rise  is  about  1000  times  the  turning  angle 

in  radians  (  I'M  =8.4).  A  10  percent  increase  in  pressure  is 

©• 

caused  by  a  0.6^^  turning  angle. 

A  comparison  is  given  in  the  following  table  between  the 
experimental  measurements  and  the  pressure  increase  predicted 
by  equations  (19)  and  (23). 
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TABLE  II 

Comparison  between  Experimental  and 
Calculated  Surface  Pressures 


Cq  X  10^ 

Ps/pood) 

Pa/Poo(2) 

0 

4.2 

.0091 

1.09 

1.08 

0 

7.6 

.0184 

1.17 

1.12 

9 

4.2 

.0191 

1.19 

1.10 

17 

4.2 

.0220 

1.23 

1.18 

25 

4.2 

.0304 

1.32 

1.30 

(1)  Calculated  from  equations  (19)  and  (23) 

(2)  Measured  data 


The  calculated  values  are  also  shown  on  figure  10,  The 
calculations  should  only  be  compared  with  the  pressures  measured 
at  the  last  four  stations,  6,  7,  8,  and  9,  because  the  experi¬ 
mental  displacement  thickness  slope  (fig.  13)  corresponds  to 
the  average  of  these  last  four  stations. 
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EXPERIMENTAL  RESULTS 

The  experimental  results  obtained  from  the  boundary  layer 
measurements  are  presented  and  discussed  in  the  following  section 
including  descriptions  of  the  Mach  number,  total  and  static  tem¬ 
perature,  and  velocity  profiles.  The  effects  of  heat  and  mass 
transfer  on  the  boundary  layer  characteristics  such  as  skin 
friction,  Stanton  number,  total  thickness,  displacement  thick¬ 
ness,  and  momentum  thickness  are  also  discussed  based  on  the  data 
contained  in  tables  III  through  XIII  which  summarize  the  data 
from  the  present  experiments. 

Mach  Number  Profiles 

A  typical  Mach  number  distribution  through  the  boundary 
layer  is  shown  in  figure  14  for  the  zero  mass  transfer  case, 
and  on  the  same  graph,  for  comparison,  are  three  mass  transfer 
cases  of  nearly  identical  Reynolds  numbers  and  wall  temperatures. 
Not  all  the  data  points  that  went  into  drawing  the  figure  are 
shown.  Zero  mass  transfer  profiles  were  obtained  for  a  relatively 
high  heat  transfer  rate  since  the  wall  temperature  is  about 
half  the  adiabatic  wall  temperature.  At  the  same  wall  tempera¬ 
ture,  the  heat  transfer  rates  are  progressively  smaller  for  the 
mass  transfer  cases. 

The  curves  siiow  a  smooth  monotonic  decrease  in  Mach  number 
from  the  free  stream  value  to  the  wall.  The  linear  region  near 
the  wall  is  roughly  about  1  mm  in  thickness  independent  of  the 
magnitude  of  the  mass  transfer  and,  hence,  also  the  boundary 
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layer  thickness.  In  addition,  the  slight  increase  in  concave¬ 
ness  vvith  mass  transfer  indicated  the  Mach  number  profiles 
cannot  be  made  similar  by  any  nondimensionalization  of  the 
distance  from  the  wall  by  the  boundary  layer  thickness. 
Temperature  Profiles 

Temperature  profiles  are  shown  in  figure  15  which  cor¬ 
responded  to  the  Mach  number  profiles  shown  in  figure  14.  Total 
temperatures,  nondimensionalized  by  the  static  temperature  at 
the  edge  of  the  boundary  layer,  were  obtained  from  equilibrium 
temperature  probe  measurements  and  show  the  strong  influence 
mass  transfer  has  on  the  total  energy  content  in  the  air  at  a 
given  distance  from  the  wall.  Although  it  is  not  very  apparent 
from  the  figure,  a  slight  increase  in  the  total  temperature 
exists  at  the  outer  edge  of  the  boundary  layer.  Conduction  of 
heat  from  relatively  hot  regions  near  the  wall  is  greater  than 
the  transport  of  energy  toward  the  wall  near  the  outer  edge, 
and  as  a  result  there  is  a  slight  accumulation  of  energy 
exceeding  that  of  the  free  stream.  High  wall  temperatures, 
i.e.,  low  heat  transfer  rates  into  the  wall,  make  the  condition 
more  pronounced. 

The  slope  of  the  total  temperature  at  each  point  of  the 
boundary  layer  is  an  indication  of  the  local  transport  of  energy, 
kinetic  as  well  as  thermal.  In  an  extended  region  (2  mm  <  y  <  6) 
the  transport  of  energy  is  increased  with  mass  transfer,  as 
evidenced  by  an  increase  in  slope  although  the  absolute  level 
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is  decreased.  Most  of  this  energy  is  absorbed  by  the  injected 
air  and  only  a  fraction  penetrates  through  the  boundary  layer 
to  the  surface.  Actually,  a  considerable  decrease  in  heat 
transfer  is  indicated  by  the  decrease  in  slope  at  the  wall 
compared  to  the  no  mass  transfer  case.  Immediately  next  to  the 
wall,  all  the  energy  is  thermal  energy  and  is  transferred  into 
the  wall  by  conduction. 

Decreasing  the  wall  temperature  at  a  given  mass  transfer 
rate  increases  the  heat  transfer  rate,  and  as  a  result  a  larger 
percentage  of  energy  is  taken  out  of  the  boundary  layer. 

Figure  16  shows  the  effect  wall  temperatures  have  on  the  tem¬ 
perature  profiles  for  zero  mass  transfer.  As  expected  the 
otal  temperature  level  in  the  boundary  layer  is  decreased  as 
the  wall  temperature  decreases,  but  the  decrease  is  not  large 
over  most  of  the  boundary  layer  (i.e.,  outer  80-90  percent). 

Most  of  the  effect  on  the  temperature  distribution  is  within 
the  10  percent  of  the  boundary  layer  nearest  the  wall. 

Static  temperature  distributions  were  calculated  from  total 
temperature  distributions  using  conventional  compressible  flow 
relations  described  previously.  All  static  temperature  curves 
approach  unity  at  the  outer  edge  of  the  boundary  layer  because 
they  have  been  nondimensionalized  by  the  free  stream  static 
temperature.  The  supply  temperature  of  the  tunnel  was  held 
constant  during  all  experiments,  and  as  a  consequence  the  static 
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temperature  in  the  free  stream  was  essentially  the  same  in  all 
cases • 

The  static  temperature  curves  have  been  connected  to  the 
value  measured  at  the  wall  oy  a  dotted  line  which  goes  through  a 
maximum  which  was  not  experimentally  observable.  This  can  be 
shown  to  be  true  because  at  the  wall  the  slope  in  temperature 
must  be  positive  for  conduction  of  heat  into  the  surface,  while 
at  the  last  measured  point  the  slope  in  static  temperature  is 
negative.  It  is  interesting  that  this  maximum  occurs  extremely 
close  to  the  surface.  Figure  17  shows  the  20  percent  of  the 
boundary  layer  nearest  to  the  surface  for  zero  mass  transfer, 
and  the  maximum  can  be  seen  to  be  about  0.1  mm  from  the  wall  or 
on  the  order  of  1  percent  of  the  boundary  layer  thickness.  In 
the  mass  transfer  cases,  the  location  of  the  maximum  is  also 
extremely  close  to  the  wall. 

Velocity  Profiles 

A  local  velocity  of  sound  can  be  calculated  from  the  static 
temperature  at  each  distance  from  the  wall  corresponding  to  a 
Mach  number  measurement  by  interpolating  between  the  measured 
temperature  data.  With  the  velocity  of  sound  and  the  Mach  number, 
the  local  velocity  can  be  computed.  Figure  18  shows  the  resulting 
velocity  distribution  through  the  boundary  layer  corresponding  to 
temperature  and  Mach  number  curves  shown  previously.  The 
velocity  is  nondimensionalized  by  the  velocity  at  the  outer  edge 
of  the  boundary  layer. 
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The  curves  have  the  typical  turbulent  boundary  layer 
appearance  with  a  sharp  decrease  in  velocity  near  the  wall, 
particularly  in  the  zero  mass  transfer  case.  The  velocity  at 
the  outer  edge  of  the  boundary  layer  is  nearly  1000  m/sec ,  and 
the  velocity  at  the  last  measured  point,  about  .3  mm  from  the 
wall,  is  420  m/sec.  Thus,  the  velocity  in  this  case  decreases 
by  40  percent  in  the  last  4  percent  of  the  boundary  layer. 

As  noted  earlier,  the  skin  friction  was  calculated  from  the 
slope  of  the  Mach  number  profile  rather  than  from  the  velocity 
profile  in  order  not  to  introduce  any  inaccuracy  from  the  tem¬ 
perature  profile  data.  However,  in  every  case  the  velocity 
slope  at  the  wall  was  computed  from  the  Mach  number  data  and  com 
pared  with  the  corresponding  velocity  curve.  Agreement  between 
the  two  procedures  was  always  within  5  or  6  percent,  which  is 
less  than  the  basic  uncertainty  in  extrapolating  either  profile 
to  zero  at  the  wall. 

Boundary  Layer  Parameters 

Based  on  the  boundary  layer  profiles  as  illustrated  by 
figures  14,  15,  and  13,  certain  properties  of  these  profiles 

Total  boundary  layer 
thickness 

Displacement  thickness 


have  been  calculated,  specifically: 

6  =  Distance  from  the  wall 

corresponding  to  u/u5  =  .99 
S 
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■=  Momentum  thickness 


ft 


-  Energy  thickness 


The  subscript  6  refers  to  the  asymptotic  value  that  any  variable 
tends  toward  near  the  outer  edge  of  the  boundary  layer.  Without 
mass  transfer,  each  of  these  parameters  increases  roughly 
linearly  with  the  distance  from  the  leading  edge  as  is  shown  in 
figure  19,  20,  21,  and  22.  The  effect  of  lowering  the  wall  tem¬ 
perature  is  to  reduce  6  and  6*  and  to  increase  6  and  ft.  The 
influence  of  6*  on  the  apparent  geometry  already  has  been 
discussed. 

Mass  transfer  also  increases  each  thickness  roughly 
linearly.  The  momentum  thickness  is  affected  the  most, 
increasing  by  a  factor  ^f  3.5  at  the  maximum  injection  rate,  and 
the  total  thickness  is  affected  the  least,  only  increasing  2.25 
times . 

The  compressible  form  factor  (displacement  thickness 
divided  by  momentum  thickness)  is  shown  in  figure  23  as  a 
function  of  blowing  rate.  These  results  show  that  heat  transfer 
has  a  significant  influence  on  the  form  factor  for  low  mass 
transfer  rates,  causing  it  to  decrease  by  almost  40  percent 
when  the  temperature  ratio  changes  from  7.6  to  4.1.  The  data 
show  that  the  form  factor  tends  to  decrease  slightly  with  mass 
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transfer  #ith  the  decrease  more  pronounced  at  the  higher  wall 
temperatures. 

Exponential  Velocity  Profile  Correlation 

One  method  of  describing  the  velocity  profiles,  which  has 
been  employed  since  the  earliest  days  of  incompressible  turbu¬ 
lent  flow  research, is  the  power-law  relation 


where  6  is  the  total  boundary  layer  thickness  and  n  is  a  para¬ 
meter  which  is,  in  general,  a  function  of  ail  the  other  parameters 
associated  with  the  flow,  i.e.,  Mach  number,  Reynolds  number, 
wall  temperature,  and  mass  transfer.  Total  thickness  data  have 
been  shown  in  figure  19  which  correspond  to  a  total  thickness 

associated  with  a  velocity  ratio  of  .99. 

Values  of  n  for  the  zero  mass  transfer  cases  are  plotted 

against  length  Reynolds  number  in  figure  24.  The  n's  were 
obtained  by  measuring  the  slope  of  the  velocity  profile  plotted 
on  log-log  paper.  A  decrease  of  n  from  between  9  and  10  to 
about  7  was  observed  with  increasing  Reynolds  number  in  the 
range  2.5  x  10®  to  10  x  10®. 

The  effect  of  mass  transfer  on  the  parameter  is  shown  in 
figure  24.  At  the  highest  injection  rate,  n  decreased  to  25 
percent  of  its  value  at  zero  injection. 

The  available  published  data  on  the  profile  exponent  were 
compiled  by  Winkler^S  and  presented  along  with  some  measurements 
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made  in  the  NOL  hypersonic  wind  tunnel.  However,  the  data  were 
presented  without  an  explanation  for  why  some  n  values  were 
higher  than  others.  The  present  data  (fig.  24)  shows  the  same 
trend  as  that  found  by  Winkler,  in  which  n  decreases  with 
increasing  Reynolds  number.  A  possible  explanation  of  these 
observations  is  that  the  high  n  values  were  obtained  in  a  transi¬ 
tional  region  even  though  the  skin  friction  and  wall  temperature 
were  a  maximum  upstream  of  the  measuring  location.  Natural 
transition  from  laminar  to  turbulent  flow  is  a  highly  unsteady 
process  and  in  the  case  of  the  present  tests,  first  breakdown 
of  the  disturbance  waves  in  the  laminar  flow  may  o  cur  15  or  more 
centimeters  ahead  of  station  6  at  one  instant,  and  in  the  next 
instant  a  laminar  boundary  layer  may  extend  downstream  of  station 
6.  The  percentage  of  time  in  which  the  boundary  layer  is  turbu¬ 
lent  increases  with  distance  until  the  turbulent  profile  is 
fully  established.  The  probe  measures  the  mean  properties  and 
cannot  detect  the  difference  between  the  normal  intermittancy  of 
the  outer  layers  of  a  turbulent  boundary  layer  and  the  momentary 
existance  of  a  laminar  boundary  layer,  the  outer  edge  of  which 
may  be  closer  to  the  wall  than  the  probe.  Therefore,  near  the 
edge  of  the  measured  boundary  layer,  the  only  difference  is  an 
increased  percentage  of  the  time  in  which  the  probe  is  in  the 
uniform  free  stream. 

Figure  24  shows  an  n  of  7  at  high  Reynolds  number,  which  is 
consistent  with  the  average  of  all  available  high  Reynolds 
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number  data  (see  fig.  10  of  ref.  35),  while  an  n  of  10  was 
observed  at  a  Reynolds  number  of  3  x  10®  in  what  otherwise 
appeared  to  be  a  turbulent  profile.  If  a  simple  averaging  is 
assumed  between  the  time  the  probe  is  in  the  free  stream  and 
the  time  in  a  n  =•  7  turbulent  boundary  layer,  then  it  can  be 
estimated  that  the  turbulent  appearing  boundary  layer  was  still 
laminar  as  much  as  30  percent  of  the  time. 


Skin  Friction  Coefficient 

The  local  skin  friction  coefficients  are  summarized  in 
tables  III  to  XIII.  Figure  25  shows  the  skin  friction  coef¬ 
ficient  plotted  against  Reynolds  number  based  on  the  distance 
from  the  leading  edge. 


where 


(28) 


Some  of  the  present  low  Reynolds  number  measurements  are  higher 
than  might  have  been  expected  for  the  fully  developed  turbulent 
boundary  layer.  Calculations  by  the  reference  temperature 
procedure^®  and  comparison  with  other  available  data37  confirm 
this  interpretation.  The  trend  at  the  higher  Reynolds  numbers, 
however,  suggests  the  fully  turbulent  condition  is  being  rapidly 


approached . 
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As  expected,  heat  transfer  increases  the  Reynolds  number 
at  the  end  of  transition.  As  can  be  seen  in  figure  25,  the 
data  at  the  highest  wall  temperature  continuously  increases  with 
decreasing  Reynolds  number  which  indicates  transition  occurs  at 
lower  values  of  Reynolds  number.  However,  at  the  lower  wall 
temperature  the  data  suggests  cf  assumes  a  maximum  value  at  the 
lowest  RCjj  which  is  characteristic  of  the  end  of  the  transition 
region.  The  maximum  of  Cf  occurs  between  station  5  and  6. 

In  figure  25,  the  Cf  data  are  plotted  against  length 
Reynolds  number  with  the  length  equal  to  the  distance  from  the 
leading  edge.  The  maximum  variation  in  Reynolds  number  is  not 
only  due  to  change  in  location  on  the  plate,  but  in  part  it  is 
due  to  a  change  in  wind  tunnel  supply  conditions.  As  a  conse¬ 
quence,  the  variation  of  transition  location  is  associated  not 
only  with  the  model  wall  temperature  but  is  also  associated  with 
the  changing  supply  pressure.  Different  supply  pressure  levels 
produce  different  free  stream  turbulence  levels,  which  in  turn 
influences  the  transition  location.  In  order  to  present  the  data 
in  such  a  way  as  to  be  independent  of  the  transition  location, 
the  Cf  data  are  plotted  as  a  function  of  momentum  thickness 
Reynolds  number  in  figure  26,  A  unique  relation  is  assumed  to 
exist  between  Cf  and  Reg  in  the  fully  developed  region,  and  the 
data  are  correlated  by  these  parameters  to  about  tne  same  extent 
as  in  the  Cf  versus  Rex  presentation. 
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Mass  transfer  decreases  the  skin  friction  by  a  considerable 
amount  as  shown  in  figure  27  where  Cf  divided  by  Cf  at  zero  mass 
transfer  (at  constant  length  Reynolds  number)  is  plotted  as  a 
function  of  Cq.  The  mass  transfer  parameter  is  twice  the  mass 
transfer  coefficient  divided  by  the  zero  mass  transfer  skin 
friction  coefficient.  The  magnitude  of  the  reduction  in  Cf  is 
emphasized  by  noting  that  an  injection  parameter  of  3.5  cor¬ 
responds  to  a  mass  transfer  rate  of  .25  percent  of  the  free 
stream  mass  flow.  Yet  this  relatively  low  mass  transfer  rate 
decreases  the  skin  friction  coefficient  to  1/4  of  its  zero 
injection  value. 

The  variation  of  momentum  thickness  with  x  provides  a 
check  on  the  consistency  of  the  skin  friction  measurements 
obtained  from  the  slope  at  the  wall  data.  The  rate  of  change  of 
momentum  thickness  with  x  can  be  related  to  Cf  through  the  inte¬ 
grated  momentum  equation.  The  relation  is. 

^  ~  ^  ^  Ca  (29) 

dX  2.  ’ 

Lines  drawn  in  figure  (28)  have  a  slope  determined  from  the 
measured  Cq  and  the  cf  obtained  from  extrapolation  of  the 
velocity  profiles  as  previously  explained.  Slopes  calculated 
in  this  way  agree  within  the  scatter  of  the  data  for  all  condi¬ 
tions  tested.  At  high  injection  rates,  the  skin  friction  plays 
only  a  small  part  in  determining  the  rate  of  change  of  0  because 
Cf/2  drops  to  only  6  percent  of  Cq. 
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Heat  Transfer  Coefficient 

At  hypersonic  speeds,  heat  transfer  is  of  greater  interest 
than  skin  friction,  and  figure  29  contains  the  results  of  the 
heat  transfer  measurements.  Both  cj  and  Stanton  number  (St) 
ratioed  to  their  respective  zero  mass  transfer  data  are  plotted 
versus  the  mass  transfer  parameter .  For  the  heat  transfer  data 
the  mass  transfer  parameter  is  formed  by  taking  the  ratio  of 
mass  transfer  coefficient  to  the  zero  mass  injection  Stanton 
number . 

Although  there  is  considerable  scatter  in  the  data,  the 
distribution  of  Stanton  number  with  mass  transfer  is  similar  to 
the  variation  of  skin  friction  coefficient.  This  leads  to  the 
conclusion  that  a  direct  correlation  exists  between  Stanton 
number  and  skin  friction  coefficient,  that  is,  a  Reynolds 
Analogy  relationship. 
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Reynolds  Analogy 

Since  both  skin  friction  and  heat  transfer  rate  were 
measured  under  the  same  conditions,  a  check  is  available  on 
Reynolds  analogy,  where  "Reynolds  Analogy"  means  in  one  case 
the  proportionality  between  Cf  and  St.  In  addition,  since  the 
velocity  and  temperature  profile  were  obtained,  the  more  funda¬ 
mental  aspects  of  Reynolds  analogy,  that  is  the  analogy  between 
momentum  and  energy  transport  at  any  point  in  the  boundary  layer, 
can  also  be  investigated. 

Figure  30  shows  the  Stanton  number  data  plotted  against 
corresponding  skin  friction  coefficients.  Zero  injection  points 
are  generally  highest  in  Cf  (filled  symbols),  and  the  mass 
transfer  points  (open  symbols)  are  the  lowest.  A  straight  line 
corresponding  to  an  ideal  Reynolds  analogy  (St  ■  Cf/2)  fits  the 
data  Quite  well,  although  the  Colburn  form  of  Reynolds  analogy 
(St  »  Cf/2Pr2/3)  is  only  a  little  higher.  The  fact  that  the 
analogy  has  been  demonstrated  in  a  hypersonic  turbulent  boundary 
layer  even  with  mass  transfer  is  a  most  significant  result. 

The  relationship  between  temperature  and  velocity  can  be 
used  to  investigate  the  local  Reynolds  analogy^  where  instead 
of  considering  the  net  flux  of  heat  from  the  boundary  layer  into 
the  wall  with  the  shear  stress  on  the  wall,  the  flux  of  energy 
is  compared  with  the  flux  of  momentum  at  each  point  in  the 
boundary  layer.  Since  local  heat  and  momentum  fluxes  are  not 
measured  in  the  boundary  layer,  the  equations  of  motion  and 
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energy  neglecting  the  pressure  gradient,  i*e., 


PU 


3X 


^  p\f  ^  -  ~ 


(29) 


PU 


_  i£ 


(30) 


are  used  to  deternine  the  coordinates  of  teeperature  and  velocity 
that  best  show  the  effect  of  local  Reynolds  analogy.  That  is. 

If  the  flux  of  energy  E  were  proportional  to  the  flow  of 
■onentuB  r,  then  from  equations  (29)  and  (30)  It  follows: 

which  is  satisfied  by  the  simple  relation: 


-  C 


(32) 


Hence,  a  linear  relationship  between  and  u  can  be  interpreted 
as  demonstrating  the  existence  of  a  local  Reynolds  analogy. 

The  constants  B  and  C  can  be  evaluated  in  various  ways.  The 
simplest  evaluations^  assumes  equation  (32)  is  valid  all  across 
the  boundary  layer.  A  sufficient  set  of  conditions  on  B  and  C 
can  be  obtained  by  using  wall  and  free  stream  conditions.  Such 
an  assumption  leads  to  the  equation: 


T,-X 


(33) 
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Equation  (33)  is  not  a  realistic  solution  of  the  energy  equation, 
but  it  does  provide  a  convenient  coordinate  systea  with  which  to 
plot  the  profile  data  in  order  to  show  the  degree  of  agreeaent 
or  disagreeaent  with  Reynolds  analogy.  If  perfect  agreeaent  is 
obtained,  then  Reynolds  analogy  exists  everywhere  in  the  boundary 
layer.  Note  that  aass  transfer  is  not  an  explicit  function  in 
equation  (29)  and  (30)  and  enters  only  through  the  boundary  con¬ 
ditions. 

Figures  31  and  32  show  the  zero  aass  transfer  profiles 

plotted  as  (Tt-Tw)/(TtrTw)  against  u/ua.  The  data  are  correlated 

0 

by  the  paraaeters  suggested  by  equation  (33) ,  but  there  is  still 
a  systeaatic  deviation  which  shows  that  a  single  linear  relation 
will  not  describe  the  aeasureaents.  In  fact,  it  is  only  reason¬ 
able  that  the  teaperature  velocity  distribution  should  bo  composed 
of  parts,  laminar  sublayer,  buffer  layer,  turbulent  layer,  etc., 
consistent  with  the  normal  picture  of  the  turbulent  velocity 
profile. 

Nevertheless,  the  temperature  velocity  profile  is  repre¬ 
sented  quite  well  by  equation  (33),  and  a  better  approximation 
is  given  by  the  empirical  equation  (34) 

=  P  ( i'l  +  ( I  -  P)  P  = 

'vv  *  *  't^  'w 

which  is  also  shown  on  figures  31  and  32, 
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Figure  33  shows  a  set  of  tenperature  velocity  profiles  at 
different  rates  of  mass  transfer.  Their  general  appearance  is 
changed  only  in  some  minor  details  by  the  mass  transfer.  For 
example,  the  zero  mass  transfer  profile  has  a  change  in  slope 
at  u/u5  -  .75  which  corresponds  to  the  location  of  the  laminar 
sublayer.  With  increasing  mass  transfer  this  point  moves  to 
lower  values  of  u/ug  consistent  with  the  lower  velocity  at  the 
sublayer  outer  edge.  However,  equation  (34)  approximates  the 
mass  transfer  temperature  velocity  relation  to  the  same  degree 
of  accuracy  as  in  the  zero  mass  transfer  cases. 
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CORRELATION  AND  ANALYSIS  OF  THE  VELOCITY  PROFILES 

Semi-empirical  analyses  of  the  turbulent  boundary  layer  are 
usually  based  on  correlations  of  velocity  profile  data.  In 
these  correlations,  the  profiles  are  somewhat  arbitrarily 
divided  into  various  layers  where  different  flow  mechanisms  are 
assumed  to  dominate.  Four  main  divisions  of  the  profile  are 
sometimes  considered,  namely,  the  laminar  sublayer,  buffer 
layer,  fully  turbulent  region  and  the  outer  wake  or  velocity 
defect  region.  These  regions  are  defined  most  simply  in  the 
nondimensional  coordinate  systems 


2.  1 

The  laminar  sublayer  was  originally  thought  of  as  a  region 
immediately  next  to  the  surface  where  viscosity  damps  out  the 
turbulent  motion.  Frequently  the  sublayer  region  is  defined  as 
extending  from  the  wall  to  a  y"^  of  5,  In  this  region  the 
velocity  gradient  is  taken  as  the  local  shear  stress  divided  by 
the  local  viscosity.  Although  the  assumption  of  a  purely 
laminar  region  has  practical  value,  it  is  not  consistent  with 
measurements^®  which  show  a  maximum  in  the  turbulence  level  at 
a  y+  of  about  5.  Nevertheless,  the  region  0  y+  5  is 
generally  referred  to  as  the  laminar  sublayer. 


N+-  ^ 
)  J 


s 


(35) 
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The  buffer  layer  is  the  transition  region  between  the 
laminar  flow  layer  near  the  wall  and  the  fully  turbulent  part 
of  the  boundary  layer.  Here  the  laminar  and  turbulent  mechanisms 
are  roughly  of  equal  Importance  in  determining  the  velocity 
distribution.  Approximately  the  region  extends  from  y"^  -  5  to 
y+  -  30. 

From  y+  -  30  to  about  y/6  -  .2,  the  velocity  profile  is 
assumed  to  be  determined  exclusively  by  the  turbulence. 

The  turbulence  does  not  extend  uniformly  to  the  ed^^e  of 
the  boundary  layer  but  becomes  intermittent  at  large  distances 
from  the  wall.  There  is  a  second  transition  region  where  the 
flow  changes  from  fully  turbulent  to  the  uniform  free  stream 
conditions.  This  region  has  the  characteristics  of  a  wake 
propagating  into  the  undisturbed  flow.  In  incompressible  flow, 
the  wake  or  velocity  defect  region  as  it  is  called  begins  at 
y/6  of  about  .2.  The  wake  region  therefore  constitutes  the 
outer  80  percent  of  the  total  boundary  layer  thickness. 

Figures  34  and  35  show  the  zero  mass  transfer  velocity 
profiles  from  the  present  experiment  in  nondimensionr 
coordinates  u+  versus  y+.  The  first  figure  contains  data  from 

the  high  heat  transfer  situation. 

Figure  35  contains  the  data  from  a  low  heat  transfer  test. 

In  both  cases  there  is  a  consistent  upward  shift  with  Reynolds 
number  in  the  fully  turbulent  region.  This  shift  is  due  to  the 
proximity  of  transition.  The  low  Reynolds  number  data  are  in  a 


t 
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region  where  fully  turbulent  conditions  are  not  yet  established. 
A  similar  conclusion  was  reached  in  an  earlier  section  in  con¬ 
nection  with  the  velocity  profile  exponent. 

As  a  consequence  of  the  transitional  nature  of  the  low 
Reynolds  number  data,  attention  will  be  exclusively  focused  on 
the  highest  Reynolds  number  tests  in  the  discussion  and  analysis 
to  follow.  The  profiles  of  main  interest  are  shown  in  figures 
36  and  37  where  velocity  profiles  for  low  ^-y^^Tjand  high 
heat  transfer  rates  are  plotted  in  u**"  versus  y+  co¬ 
ordinates.  Although  the  general  level  of  all  the  profiles  is 
the  same,  there  is  a  definite  Increase  in  slope,  which  may  be 
attributed  to  the  increased  mass  transfer. 

M ixing  Length  Theory  of  the  Fully  Turbulent  Region 

The  fully  turbulent  region  is  the  most  significant  part  of 
the  boundary  layer  and  determines  much  of  its  appearance,  even 
if  it  does  not  always  constitute  a  physically  large  part  of  the 
profile.  Prandtl  proposed  that  the  turbulent  shear  stress 
could  be  written  as  follows: 


Essentially,  the  formulation  is  based  on  dimensional  reasoning 
and  a  very  simplified  physical  model  of  turbulence.  Furthermore, 
since  the  fully  turbulent  layer  occurs  close  to  the  wall,  it 
was  assumed  that  the  shear  stress  was  constant  and  for  all 
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practical  purposes  equal  to  the  wall  value.  The  value  of  the 
mixing  length.  I  ,  was  evaluated  not  fro*  theory  but  was  based 
on  observation  of  experimental  data.  The  mixing  length,  it  was 
suggested,  was  proportional  to  the  distance  from  the  wall,  and 
the  constant  of  proportionality  was  found  to  be  a  universal 
constant . 

In  incompressible  flow,  these  relations  can  be  manipulated  into 
the  "log-law"  velocity  profile,  the  slope  of  which  is  given  by 

In  -NlikZf.  (36) 

d  y  ky 


(37) 


or  in  nondimensional  coordinates: 

d_u^  _ 

a 

Equation  (37)  can  be  Integrated  to  give  the  log-law  which  Is 

-i-C 

The  constants  k  and  C  have  been  determined  from  experiments,  as 


k  -  .4  and  C  “  5,5+  .5,  for  a  smooth  flat  plate. 

Effects  of  compressibility  are  generally  accounted  for  in 
the  shear  stress  equation  (1)  by  using  the  local  density. 
Effects  of  Mass  Transfer  on  Fully  Turbulent  Region 

Mass  transfer  into  the  boundary  layer  changes  the  velocity 
gradient  relations  by  adding  a  terra  to  the  local  shear  stress. 
The  shear  stress  is  not  a  constant  but  increases  with  velocity 
and  mass  transfer  as  can  be  shown  by  consideration  of  the  moraen 
turn  and  continuity  equations. 


(15) 


Momentum  Equation 

Continuity  Equation 

A  (15) 

3x 

It  has  been  shown  by  experiment  that  when  the  boundary  layer 
is  thin,  the  first  term  in  equations  (29)  and  (15)  which  contain 
derivatives  with  respect  to  x  can  be  neglected  when  compared  to 
the  other  terms.  (See  discussion  on  page  6.)  The  continuity 
equation  (15)  can  then  be  integrated  immediately  to  give  plTs 
constant.  The  constant  is  evaluated  from  the  known  condition  at 
the  wall  where  the  mass  flow  is 

If  the  first  term  in  equation  (29)  is  neglected  and  the 
pV  in  the  second  terra  is  set  equal  to  then  it  can  also 

be  integrated  from  the  wall  (y  “  0,  u  «  0,  arbi¬ 

trary  point  (y,  u.T). 

(pu')^u  =  r-r^  (39) 


(39) 


or 

t  =  (40) 

Using  equation  (40),  Prandtl’s  equation  for  the  slope  of  the 
turbulent  velocity  profile  including  the  effect  of  mass  transfer 


becomes ; 


^  |fw/P 


c<l 


(41) 
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Transformation  of  Compressible  Velocity  Profile  (with  Heat  and 
Mass  Transfer)  into  the  Incompressible  Velocity  Profile 

The  density  ratio  in  equation  (41)  can  be  expressed  in 
terms  of  the  local  velocity  through  the  Crocco  temperature- 
velocity  relation,  and,  thus,  it  is  convenient  to  transfer  all 
velocity  terms  to  the  left  hand  side  of  the  equation.  A  trans¬ 
formed  velocity  coordinate  can  then  be  introduced  which  satisfies 
the  following  equation: 


\ 


u 


(42) 


Equation  (41)  in  terms  of  u’*’’^  and  y+  has  the  same  form  as  the 
incompressible  zero  mass  transfer  velocity  slope  and  has  the 
same  solution. 


U++=  +c 


f  A’i\ 


Equation  (43)  applies  only  in  the  fully  turbulent  part  of  the 
profile  and  the  constant,  C,  must  be  obtained  from  experiments. 
Results  of  the  compressible  but  zero  mass  transfer  analysis  of 
Van  Driest^,  Harkness^®,  and  Moore^^  can  be  produced  from  the 
above  transformation.  In  each  of  these  analyses,  a  density- 
velocity  relationship  was  introduced,  and  equation  (42)  or  its 
equivalent  was  integrated  to  obtain  a  complex  expression  in 
terms  of  the  velocity  for  the  left  hand  side  of  equation  (43). 
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For  example,  Van  Driest  used  the  following  form  for  the  density 


ratio: 


where 


(44) 


Equation  (42)  can  be  integrated  in  this  case  to  give 


A 


e. 


=  ii«v|+  +  C 


(45) 


Although  the  right  hand  side  implies  the  equation  is  valid  only 
in  the  fully  turbulent  region,  the  left  hand  side  has  been  inte¬ 
grated  from  the  wall  (u/uj  -  0)  to  some  arbitrary  point  in  the 
turbulent  region. 

Another  example  is  the  result  obtained  by  Moore^^  which 
was  originally  derived  starting  from  von  Karman's  equation. 
However,  the  same  result  can  be  arrived  at  through  use  of 
equation  (42).  Moore  employed  the  following  relation  for  the 
density  ratio,  which  corresponds  to  equation  (44)  when  the  wall 
temperature  equals  the  adiabatic  wall  temperature 


(46) 
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Equation  (42)  and  (43)  become 


a 


-i 


\fh^- 

T.  U 


w  ‘*8 


-  i-  C 

Vc 


(47) 


The  more  general  case  including  mass  injection  requires  the 
integration  of  equation  (42)  combined  with  the  density  ratio 
given  by  equation  (44).  Such  an  integration  has  been  performed 
by  Rubesin^^  where  equation  (42)  is  put  in  the  following  form 

/U/Mg 


A  “ 


(48) 


The  integration  gives 


a 


where 
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£■  =  ( 6  -  'Ib^+4-a^  )/z 
p=  (b- 
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Elliptical  Integral  of  first 


kind 


(a)  if 
then 


^  -1 


I _ BL 


Ti-  S 


NOLTR  64-99 


(b)  if 


then 


<1 

u 


b-E 


As  pointed  out  above,  the  integration  of  equations  (42)  and  (43) 
is  restricted  to  the  fully  turbulent  part  of  the  boundary  layer. 
Nevertheless,  the  equations  cited  of  Van  Driest  and  Moore  were 
integrated  from  the  wall,  a  point  not  in  the  turbulent  part  of 
the  profile.  Moore  has  even  suggested  that  the  right  hand  side 
of  equation  (43)  be  replaced  by  a  law  of  the  wall  (e.g..  Coles, 
ref.  39)  with  a  term  added  to  account  for  the  variation  of  C 
with  heat  transfer.  Since  the  added  term  cannot  be  Just  a 
function  of  the  heat  transfer  rate  and  have  the  equations  cor¬ 
rectly  transformed  in  the  sublayer  and  at  tht  wall,  it  appears 
that  such  a  procedure  is  not  completely  sufficient  to  reduce  the 
compressible  profile  with  heat  transfer  to  the  incompressible 
profile.  In  the  data  reduction  for  the  following  analysis,  it 
will  be  assumed  that  the  velocity  coordinate  u’^ 


2  ^ 


(50) 


is  a  function  of  the  y**"  coordinate,  and  it  will  be  the  objective 


to  determine  from  the  available  data  the  variation  of  the 
profile  parameter  C  with  heat  transfer,  Mach  number,  and  mass 
transfer . 
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Prof ile  Constant  Determination 

There  is  insufficient  understanding  of  the  turbulent 
boundary  layer  to  be  able  to  calculate  the  profile  constant  in 
the  "log- law"  from  theory.  The  value  of  the  constant,  C,  may  be 
considered  as  determined  by  the  lower  limit  (of  velocity  and 
distance  from  the  wall)  of  the  fully  turbulent  region,  or  it 
may  be  considered  as  the  upper  boundary  of  the  laminar  sublayer- 
buffer  layer  regions.  Near  the  wall,  the  velocity  is  determined 
by  the  distribution  of  laminar  and  turbulent  shear  stresses, 
and,  therefore,  the  constant  C  is  also  determined  by  the  distri¬ 
bution  of  these  shear  stresses.  Some  analyses  have  simplified 
this  complex  problem  by  neglecting  the  buffer  region  and  only 
considering  a  laminar  sublayer  and  a  fully  turbulent  region 
where  the  flow  changes  from  laminar  to  turbulent  discontinuously , 
The  problem  of  determining  the  location  of  the  laminar-turbulent 
interface  still  remains,  and  various  assumptions  have  been  tried 
in  the  analyses.  However,  the  interface  location  is  not  a 
physically  real  point  in  the  flow,  and,  therefore,  it  can  be 
questioned  whether  anything  fundamental  can  be  obtained  by  cor¬ 
relating  data  using  that  concept.  The  constant,  C,  on  the 
other  hand,  provides  a  convenient  way  of  correlating  the  data, 
and  it  may  be  thought  of  as  referring  to  a  physical  location  on 
the  boundary  layer  profiles. 
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Incompressible  Flow 

The  first  measurements  of  the  profile  constant  were 
obtained  in  incompressible  pipe  flow^®,  and  over  the  years  a 
large  amount  of  incompressible  data  has  been  accumulated. 

Despite  the  number  of  measurements,  there  is  still  a  rather 
large  uncertainty  in  the  reported  values.  For  example,  Coles^S 
in  a  review  of  selected  incompressible  flat  plate  boundary  layer 
measurements  has  concluded  that  C  has  a  value  of  5.10+  .2  at 
Reynolds  numbers  above  2000,  Ross^^  investigated  a  large  number 
of  sources  in  which  data  were  reported  for  various  configura¬ 
tions  (pipes,  channels,  flat  plates,  and  airfoils)  and  found 
C  -  5.55+  .4.  Using  much  of  the  same  data,  Clauser'*^  obtained 
C  -  4.9.  The  fact  that  C  has  not  been  measured  more  precisely 
during  the  past  30  years  gives  some  indication  of  the  difficulty 
of  performing  these  experiments.  Even  recent  data  show  widely 
spread  results.  For  example.  Smith  and  Walker^^  determined  a 
value  of  C  -  7.15  fitted  their  data  best  at  high  Reynolds 
numbers,  whereas  Reynolds  et  al"^"*  obtained  C  «  4.4. 

Compressible  Adiabatic  Flow 

Coles^^  has  shown  that  when  the  adiabatic  compressible 
velocity  profiles  are  transformed  into  coordinates  similar  to 
those  presented  here,  they  reduce  approximately  to  the  incorapres 
sible  profiles.  This  is  equivalent  to  saying  the  constant,  C, 
has  about  the  same  value  as  in  low  speed  flow.  Figure  38  shows 
the  few  data  points  available  for  the  adiabatic  supersonic 
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case37, 45-48  as  well  as  a  few  measurements  of  the  incompressible 
boundary  Most  of  this  data  falls  between  4.5 

and  6.5,  and  they  are  apparently  independent  of  Mach  number 
below  M  -  5.  The  only  data  at  higher  Mach  numbers  are  those  of 
Matting,  et  al^^,  which  were  obtained  in  helium.  However, 
because  of  the  uncertainty  in  Interpreting  the  results  in  helium 
and  because  of  Inconsistencies  in  those  particular  profiles, 
Matting's  data  were  not  included  here. 

Compressible  Flow  with  Heat  Transfer 

Harkness^®  and  Moore^^  have  considered  the  effect  of  heat 
transfer  on  the  compressible  turbulent  boundary  layer.  Although 
they  employed  the  von  Karman  equation  and  presented  the  results 
in  terms  of  the  sublayer  thickness,  their  conclusions  apply 
equally  well  to  the  value  of  C.  They  determine  the  effect  of 
heat  transfer  from  the  experimental  data  of  Monaghan  and  Cooke^®, 
Lobb,  Winkler  and  Persh47,  ^nd  Winkler  and  Cha35.  From  these 
experiments,  it  was  found  that  C  increases  quite  rapidly  with 
heat  transfer  into  the  wall.  The  data  used  by  Harkness  and 
Moore  has  been  reanalysized  to  obtain  the  profile  constant  and 
figure  39  shows  C  as  a  function  of  the  nondimensional  wall  tem¬ 
perature,  which  is,  of  course,  a  direct  indication  of  the  heat 
transfer. 

The  highest  Reynolds  number  profiles  of  the  present  experi¬ 
ments  have  been  transformed  through  the  use  of  equation  (42), 
and  the  no  blowing  measurements  are  shown  in  figure  40.  They 
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show  the  expected  trend  of  increasing  C  with  increasing  heat 
transfer.  The  values  of  C  obtained  from  these  curves  are  also 
plotted  in  figure  39  and  are  generally  consistent  with  the 
earlier  results.  It  appears  from  the  figure  that  the  high  heat 
transfer  profile  constant  is  a  little  high  but  still  within  the 
scatter  of  the  available  data.  Since  the  data  support  an 
increase  of  C  with  decreasing  wall  temperature,  it  is  interesting 
to  speculate  on  the  more  fundamental  reasons  for  such  behavior. 

As  indicated  earlier,  C  is  determined  by  the  distribution  of 
laminar  and  turbulent  shear  stresses  in  the  sublayer  and  buffer 
layer.  If  the  contribution  of  the  turbulent  shear  were  decreased, 
then  laminar  conditions  would  contribute  more  in  determining 
the  local  velocity.  Since  the  effect  of  turbulence  is  to 
decrease  the  velocity  gradient  from  the  higher  laminar  velocity 
gradient,  it  follows  that  a  decrease  in  turbulence  would  result 
in  a  higher  velocity  gradient.  Since  the  measured  velocity  at 
the  lower  edge  of  the  fully  turbulent  region  increased  under 
conditions  of  heat  transfer  into  the  wall  compared  to  the  adia¬ 
batic  case,  it  appears  that  removing  energy  from  the  boundary 
layer  lowers  the  turbulence  level  near  the  wall,  and  conversely 
heat  transfer  away  from  the  wall  tends  to  increase  the  turbulence 
level . 

As  is  well  known,  heat  transfer  produces  just  such  an  effect 
on  the  stability  and  transition  of  a  laminar  boundary  layer^®. 
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Thus,  there  appears  to  be  an  analogous  effect  in  fully  turbulent 
flow. 


It  is  not  clear  whether  there  is  a  change  in  the  distribu¬ 
tion  of  turbulence  or  whether  there  is  an  increased  thickness  of 
a  purely  laminar  region  or  both.  For  example,  if  a  purely 
laminar  region  exists,  a  change  in  its  thickness  could  be 
reflected  directly  as  a  change  in  the  profile  constant.  A 
stability  criterion  might  be  found  for  such  a  case,  which 
reflects  the  effect  of  heat  transfer,  similar  to  the  criterion 
suggested  by  Van  Driest  for  the  location  of  transition  of  the 
laminar  turbulent  interface^^. 

Although  a  reasonable  correlation  of  the  data  on  the  basis 
of  the  temperature  ratio  is  shown  in  figure  39,  the  following 
correlating  parameter  was  also  considered: 


4-^  \  Te 
^  2-  Te-T^ 


(51) 


Equation  (51)  is  derived  from  the  Crocco  temperature-velocity 
relation  by  hypothesizing  that  the  change  in  sublayer  thickness 
is  related  to  the  thickness  of  the  region  between  th'  ill  and 
the  maximum  in  the  static  temperature.  As  was  pointed  out  when 
the  temperature  profiles  were  discussed,  the  maximum  static 
temperature  occurs  very  near  the  surface  and  in  a  region  where 
the  nondiraensional  distance  y'*’  is  approximately  equal  to  u+. 
Thus,  equation  (44)  can  be  differentiated  and  set  equal  to  zero 
to  find  the  0+  -  y+  point  corresponding  to  the  maximum 
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temperature.  The  result  is  given  by  equation  (51).  The  varia¬ 
tion  of  C  as  a  function  of  A+  Is  shown  in  figure  41.  The  cor¬ 
relation  is  slightly  better  than  that  obtained  in  figure  39. 

The  Effect  of  Mass  Transfer  on  C 

Two  analyses  may  be  cited  which  used  a  similar  approach  to 

that  used  here  for  describing  the  effect  of  mass  transfer  on 
the  boundary  layer.  Both  the  work  of  Rubesinl^  and  Dorrance 
and  DorelS  start  with  the  same  equations  and  derive  a  relation 
which  essentially  corresponds  to  equations  (42)  and  (43),  At 
the  time  of  their  analyses  (1954),  no  velocity  profile  data  were 
available  for  mass  transfer  into  the  compressible  turbulent 
boundary  layer.  Dorrance  and  Dore  overcame  this  difficulty  by 
making  various  mathematical  assumptions  in  deriving  a  skin 
friction  relation.  They  then  evaluated  the  constants  k  and  C 
from  the  incompressible,  zero  mass  transfer  limiting  case  using 
the  Karman-Schoenherr  equation.  From  these  and  some  additional 
assumptions,  they  arrived  at  a  value  of  k  =  .393,  and  C  was 
taken  as  either  4.13  or  6,53, 

In  effect,  Rubesin's  approach  was  to  determine  the  constant 
C,  by  defining  a  sublayer  thickness,  where  u+a  “  y'‘'a» 
inner  edge  of  the  fully  turbulent  part  of  the  profile.  He  used 
an  equation  similar  to  (43)  to  describe  tbs  velocity  between  the 
sublayer  and  the  free  stream.  A  further  assumption  was  needed 
to  define  the  sublayer  thickness,  and  it  was  assumed  that  incom¬ 
pressible  data  without  heat  or  mass  transfer  could  be  used  to 
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determine  the  point  y^a' * 

k  -  .392  were  recommended  which  corresponds  to  a  C  of  6.5. 

The  present  experimental  data  provide  an  opportunity  to 
extend  the  experimental  data  on  the  constant  C  to  include  the 
effect  of  mass  transfer.  Figures  42  and  43  show  the  result  of 
transforming  the  velocity  coordinate  according  to  equation  (42) 
and  using  the  experimental  temperature  profile  and  wall  mass 
transfer  rate.  Both  sets  of  profiles  at  two  wall  temperatures 
show  the  same  trend.  As  the  mass  transfer  rate  is  increased, 
the  profiles  are  lowered,  and,  as  a  consequence,  C  also  decreases. 

It  is  very  interesting  to  note  that  a  semi-log  relationship  fits 
the  measured  points  very  well  in  the  fully  turbulent  part  of 
the  boundary  layer  from  y+  -  30  to  100.  There  is  no  difficulty 
in  drawing  a  line  of  slope.  1/k.  through  the  points  and  determining 
C.  Figure  44  summarizes  the  results  by  showing  the  change  in  C 
with  the  dimensionless  mass  transfer  rate.  c<j.  Both  wall  tem¬ 
perature  conditions  approximately  exhibit  the  same  decrease  of 

C  with  increasing  injection  coefficient. 

The  only  other  experimental  data  on  the  turbulent  velocity 
profiles  with  mass  transfer  are  the  measurements  reported  by 
Mickley  and  DavisSl  and  Tewfik^S  incompressible  flat  plate 

boundary  layers.  Mickley’ s  measurements  Indicate  the  opposite 
trend  from  the  present  M  -  6.7  data  and  shows  a  strong  Reynolds 
number  effect.  Tewfik’s  data,  on  the  other  hand,  shows  a 
marked  decrease  in  C  for  the  high  injection  rate  profile  given 
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in  reference  52.  Obviously,  more  profile  measurements  are 
required  to  establish  the  variation  of  C  with  Mach  number  and 

possibly  the  effect  of  surface  material. 

Theoretical  analyses  of  mass  transfer  consider  the  idealized 
case  of  perfectly  uniform  injection.  In  practice,  uniform 
Injection  can  only  be  approximated  with  semi-solid  materials, 
and,  therefore,  it  is  not  surprising  that  mass  transfer  has  a 
destabilizing  effect  on  the  sublayer.  Although  the  injection 
of  air  has  been  shown  to  reduce  the  stability  of  laminar  boundary 
layers53,  such  a  phenomenon  would  not  explain  the  large  decrease 
in  C.  If  just  the  stability  of  the  laminar  sublayer  were 
Involved,  the  decrease  in  C  would  be  expected  to  asymptotically 
approach  some  value  that  would  correspond  to  the  fully  turbu¬ 
lent  conditions  throughout  the  sublayer.  It  appears  more  likely 
that  the  destabilizing  effect  is  due  to  local  nonuniformity  of 
the  injection,  and  the  effect  is  similar  to  that  of  roughness. 
Although  the  profile  constant  is  given  by  Schlichting^^  in  a 
somewhat  different  form,  examination  shows  that  C  decreases 
with  increasing  roughness  height  in  incompressible  flow.  For 
small  roughness  heights,  the  decrease  in  C  is  explained  by  an 
increased  turbulence  level  in  the  region  near  the  wall  which  is 
related  to  the  vortices  shed  from  the  disturbance  elements. 

The  present  porous  plate  is  composed  of  small  spheres  40 
microns  and  smaller  in  diameter  compressed  under  high  pressure 
until  a  flat,  smooth  surface  was  obtained.  The  flow  passages 
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through  the  material  results  in  the  air  being  injected  as  tiny 
jets  of  probably  unequal  strength  and  almost  random  direction. 

From  a  macroscopic  view  point,  the  mass  transfer  is  very  uni¬ 
form,  but  from  a  more  microscopic  view  point,  the  jets  may  be 
expected  to  induce  a  higher  level  of  turbulence  in  the  sublayer 

and  buffer  layer. 

Skin  Friction  Law 

Empirically  correlated  velocity  profiles,  such  as  the 
”log-law,"  are  also  essentially  skin  friction  laws  because  they 
contain  the  skin  friction  as  a  parameter.  The  conversion  to  a 
skin  friction  Reynolds  number  relation  is  obtained  by  inserting 
the  velocity  and  temperature  profile  correlations  into  the  inte¬ 
grated  momentum  equation.  The  skin  friction  law  is,  of  course, 
only  as  accurate  as  the  accuracy  of  the  correlations.  On  the 
other  hand,  because  it  is  an  integral  method,  the  result  does 
not  require  an  exact  solution  near  the  wall.  Correlated  veloc¬ 
ity  and  temperature  profiles  developed  in  the  preceeding  analysis 
were  programmed  on  a  high  speed  digital  computer,  and  the  skin 
friction  was  computed  from  the  integral  momentum  equation  as  a 
function  of  the  momentum  thickness  Reynolds  number.  The  general 
availability  of  computing  machines  eliminates  the  need  for 
closed  form  solutions  to  the  equations  which  normally  involve 
some  mathematical  approximations.  The  equations  which  were 
programmed  for  the  U.  S.  Naval  Ordnance  Laboratory  IBM  7090 

computer  were: 
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Momentum  thickness  (definition) 

/  6 


(52) 


Velocity  profile  ("log-law”) 

.w•/U^  ,, 


(53) 


Q+-L  yi+  ^  _  L  jU\j^  -V  c 

^  Ti  2^  U  ^ 

Uj 

Density  profile  (modified  Crocco  temperature- velocity 


(42) 


equation) 


'i;  -  X  -  \  +  -'kX*  X  _  Te-'T* 


u\2- 


“s 


(44) 


The  parameters  specified  for  each  calculation  were: 

Cr  C  = 

f  ^  ' '  Ta  '  =r  '  ^ 

—  was  adjusted  to  account  for  mass  transfer  effects  using 
figure  8.  Heat  transfer  and  mass  transfer  effects  on  the  profile 
constant,  C.  were  calculated  using  the  following  approximate 
expression,  which  is  a  reasonably  good  fit  to  the  data  shown  in 


figures  41  and  44, 


(54) 


74 


NOLTR  G4-99 


Although  equation  (42)  can  be  integrated  in  a  general  way  in 
terms  of  elliptical  integrals,  it  is  simpler  to  perform  the 
operation  by  numerical  machine  methods.  In  the  calculation, 
only  the  "log-law"  was  considered,  and  the  sublayer,  buffer 
layer,  and  wake  region  were  neglected.  Thus,  the  momentum 
thickness  Reynolds  number  may  be  somewhat  higher  than  would  be 
the  case  if  these  regions  were  included.  However,  no  mathemat¬ 
ical  approximation  is  required  to  obtain  the  momentum  thickness. 

The  velocity  ratio  was  considered  the  independent  variable, 
and  at  fixed  values  of  U/U5,  calculations  were  made  of  the  tem¬ 
perature  ratio,  y'*’,  and  the  value  of  the  momentum  thickness 
integral  out  to  that  0/05  point.  Adequate  accuracy  was  obtained 

using  increments  in  u/u^  of  ,04. 

Two  sets  of  calculations  corresponding  to  two  wall  tem¬ 
perature  ratios  of  7  and  4  are  shown  in  figures  45  and  46,  For 
each  wall  temperature,  the  zero  mass  transfer  and  three  mass 
transfer  cases  are  shown.  Relatively  good  agreement  is  obtained 
between  the  calculations  and  the  experimental  data.  Lowering 
tht  wall  temperature  ratio  from  7  to  4  corresponds  to  a  change 
in  wall  temperature  of  about  1S0®C  and  results  in  a  decrease  in 
the  skin  friction  coefficient  of  11  percent.  It  is  a  relatively 
simple  addition  calculation  to  find  a  length  Reynolds  number 
which  corresponds  to  the  momentum  thickness  Reynolds  number. 

The  integral  momentum  equation  can  be  put  in  the  following  form: 
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However,  the  main  problem  is  in  interpreting  the  length,  x,  in 
terms  of  the  physical  length  from  the  leading  edge.  An 
effective  starting  length  for  a  fully  turbulent  boundary  layer 
must  be  defined,  and  since  the  objective  of  showing  the  cor¬ 
respondence  between  the  velocity  profile  correlation  and  skin 
friction  law  has  been  accomplished,  further  consideration  of 
the  skin  friction  coefficient  length  relation  Reynolds  number 
should  await  a  better  understanding  of  transition  from  laminar 
to  turbulent  boundary  layers. 
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CONCLUSIONS 

Skin  friction,  heat  transfer,  and  boundary  layer  character¬ 
istics  have  been  measured  on  a  flat  plate  at  a  Mach  number  of 
6.7,  Reynolds  numbers  between  3  x  10®  and  11  x  10®,  wall  to 
free-stream  temperature  ratios  between  4  and  8,  and  with  various 
rates  of  air  transpiration.  The  experimental  boundary  layer 
profiles  are  presented  and  analyzed  from  the  viewpoint  of  mixing 
length  theory.  The  main  conclusions  are  as  follows: 

1.  Skin  friction  and  heat  transfer  rates  have  been  found  to 
decrease  with  air  transpiration  at  Mach  number  6.7  in  about  the 
same  amount  as  has  been  found  by  other  investigators  at  lower 
Mach  numbers.  An  injection  rate  of  one  quarter  of  one  percent 
of  the  free-strea*  mass  flow  caused  the  skin  friction  to  be 
r^educed  to  twenty-five  percent  of  its  zero  injection  value. 

2.  The  total  boundary  layer  thickness  and  the  displacement 
thickness  decreased  when  the  wall  temperature  was  lowered,  but 
the  momentum  and  energy  thicknesses  increased.  Air  injection 
into  the  boundary  layer  was  accompanied  by  increases  in  all  of 
these  thicknesses,  although  the  momentum  and  energy  thicknesses 
increased  more  rapidly  than  ths  total  boundary  layer  thickness. 

3.  Both  skin  friction  and  heat  transfer  rates  were  measured 
under  essentially  identical  conditions,  and  thus  it  was  possible 
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to  experimentally  verify  Reynolds  analogy.  Stanton  numbers  were 
found  to  be  approximately  equal  to  one-half  of  the  corresponding 
skin  friction  coefficients.  This  was  found  valid  both  without 
and  with  air  transpiration  into  the  boundary  layer.  Colburn's 
form  of  Reynolds  analogy  predicts  that  the  Stanton  number  is 
equal  to  about  63  percent  (S^  -  cf/2Pr2/3,  Pr  =  .7)  of  the  skin 
friction  coefficient  which  is  somewhat  higher  than  the  percentage 
found  in  the  present  results  but  still  within  the  scatter  of  the 
data. 


4.  Temperature  and  velocity  profiles  were  also  obtained 
under  essentially  identical  conditions,  stnd  these  profile 
mfisureroents  show  that  the  modified  Grocco  equation  adequately 
describes  the  relationship  between  total  temperature  and  velocity 
within  the  turbulent  boundary  layer.  This  conclusion  was  found 
to  apply  to  boundary  layers  without  and  with  air  transpiration. 
The  modified  Crocco  equation  can  be  written  as  follows: 


Temperature  profiles  calculated  from  the  above  equation  do  not 
describe  all  the  features  of  the  measured  temperature  profile, 
and,  therefore,  it  should  only  be  considered  as  a  good  approxi 
mation. 
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5.  Prandtl's  mixing  length  analysis  has  been  shown  to 
acceptably  describe  the  experimental  velocity  profile  in  the 
fully  turbulent  region  of  the  boundary  layer  if  the  effects  of 
heat  and  mass  transfer  on  the  laminar  sublayer  are  taken  into 
account.  The  mixing  length  theory  was  used  to  develop  the  fol¬ 


lowing  transformation  for  the  local  velocity. 


= 


4  i 


I  .  o  U.* 

I  +  2  Uj 


(42) 


This  transformed  velocity  is  related  to  the  distance  from  the 
wall  by  the  same  semi- logarithmic  equation  established  for 


incompressible  turbulent  boundary  layers. 

IA'*'”*’  =  *1^  Jl/Yxy’*'  -V  C  (43) 

\  Unlike  incompressible  flow,  however,  the  constant  C  is  found 
f^om  the  present  and  other  available  data  to  be  a  function  of 
hfeat  and  mass  transfer,  but  based  on  the  available  data,  the 
constant  is  independent  of  Mach  number  under  adiabatic  wall  con¬ 
ditions  and  equal  to  its  incompressible  value  of  5.5. 


6.  Heat  transfer  was  found  to  increase  or  decrease  C 
depending  on  whether  the  heat  transfer  is  toward  or  away  from 
the  wall.  Apparently,  the  change  in  C  with  heat  transfer  is 
associated  with  a  change  in  the  turbulence  level  in  the  laminar 
sublayer  which  is  similar  to  the  effect  of  heat  transfer  on 
laminar  boundary  layer  stibility.  An  increasing  rate  of  heat 
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transfer  into  the  wall  is  accompanied  by  an  increasing  laminar 
sublayer  thickness.  This  fact  suggested  that  the  distance 
between  the  wall  and  the  maximum  in  the  static  temperature 
profile  might  be  a  suitable  length  with  which  to  correlate  the 
change  in  C.  The  maximum  in  the  temperature  profile  occurs  very 
near  the  surface  inside  the  laminar  sublayer,  and  it  is  directly 
related  to  the  heat  transfer  rate.  The  present  and  other 
published  profile  constant  data  were  correlated  rather  well  by 
the  nondimensional  distance  to  the  temperature  peak, 

.+ 


u 


(51) 


which  was  calculated  for  each  case  using  the  modified  Crocco 
temperature  velocity  relation. 


7.  Air  transpiration  into  the  turbulent  boundary  layer 
decreases  the  profile  constant  for  the  conditions  tested  in  a 
manner  analogous  to  roughness.  It  should  be  noted,  however, 
that  a  large  part  of  the  effect  of  mass  transfer  on  the  velocity 
profile  is  taken  out  by  the  transformation,  and  the  change  in 
the  profile  constant  is  an  additional  effect.  Apparently, 
transpiration  increases  the  turbulence  level  in  the  laminar 
sublayer  and  decreases  the  sublayer  thickness.  This  may  be 
related  to  the  nonuniformity  of  the  mass  transfer  process  when 
viewed  on  the  scale  of  the  pores  in  the  porous  plate  or  the 
thickness  of  the  sublayer.  It  was  found  that  the  effects  of 
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heat  and  nass  transfer  on  the  velocity  profile  constant  could 
be  accounted  for  by  the  following  approxinate  equation: 


C  =5.5  +  .'2.7S^^ut  +-.59M0^c<^ 


(54) 


8.  Skin  friction  coefficients  were  coeputed  as  a  function 
of  BoeentuH  thickness  Reynolds  nunber  by  nueerical  techniques 
using; 

a.  The  definition  of  the  Bomentun  thickness 

b.  The  transforeed  velocity  profile  ("log- 
law*'  neglecting  the  sublayer »  buffer 
layer,  and  wake  region) 

c.  The  Bodlfied  Crocco  tewperature-veiocity  relation 

d.  An  eepirlcal  profile  constant  accounting  for 
the  effect  o2  heat  and  mass  transfer  derived 
froB  the  experisental  data 

Good  agreenent  was  obtained  between  the  Beasurenents  this 

siBple  analysis. 

In  addition  to  the  principal  conclusions,  the  following 
points  appeared  to  be  significant  particularly  for  those 
interested  in  turbulent  boundary  layer  experimentation  at 
hypersonic  Mach  nuBbers. 
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9.  At  hypersonic  speeds,  boundary  layer  flow  field  inter¬ 
action  Is  iaportant  in  the  calculation  of  the  boundary  layer 
characteristics  especially  when  mass  transfer  is  involved.  The 
pressure  rise  in  the  boundary  layer  was  found  to  be  given 
approxieately  by 


where Ca+‘*-*repre8eut8  the  angle  through  which  the  flow  nust  turn 
“  dv 

to  accommodate  to  the  growing  boundary  layer. 

10.  The  experiments  described  here  were  measured  in  a 
turbulent  boundary  layer  behind  natural  transition.  However, 
the  approach  to  the  fully  turbulent  condition  was  found  to  be 
very  slow,  and  it  was  necessary  to  go  to  the  highest  possible 
Reynolds  number  per  foot  in  the  NOL  Hypersonic  Tunnel  No,  4. 

The  observed  slow  approach  to  fully  developed  turbulent  condi¬ 
tions  is  believed  to  be  related  to  the  unsteady  movement  of 
transition  which  at  a  given  time  may  be  a  relatively  largo 
distance  up  or  downstream  of  the  mean  transition  location.  As 
a  consequence,  some  measuring  stations  that  were  initially 
thought  to  be  in  the  fully  turbulent  region  exhibited  fuller 
velocity  profiles  (i.e.,  lower  exponents  in  the  power  law 
velocity  correlation  than  expected  for  fully  turbulent  condi¬ 
tions).  More  research  nsucds  to  be  done  to  determine  the 
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characteristics  of  this  part  of  the  transition  region  and 
on  developing  techniques  for  tripping  the  boundary  layer  at 

hypersonic  speeds. 
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APPENDIX  I 

ACCURACY  OF  PITOT  PROBE  MEASUREMENTS 

Recently  the  use  of  iapact  pressure  probes  has  been  criti¬ 
cized  because  of  several  possible  sources  of  inaccuracies 
associated  with  the  instrueent.  Evidence  of  a  contradiction 
between  skin  friction  aeasureaents  obtained  by  the  skin  friction 
balance  technique  and  Pitot  pressure  surveys  has  been  reported 
by  Matting,  et  al^^,  for  hypersonic  turbulent  boundary  layers. 

No  satisfactory  explanation  for  the  discrepancy  has  yet  been 
offered.  Therefore,  the  several  possible  sources  of  inaccuracy 
are  considered  here  with  the  objective  of  estisating  under  what 
conditions  and  to  what  extent  the  probe  gives  erroneous  results. 
Probe  Geosetry— First  of  all,  a  round  Pitot  tube  was  chosen 
because  of  the  greater  certainty  in  establishing  the  position 
of  the  probe  with  respect  to  the  wall.  An  uncertainty  exists 
when  two-disenslonal  probes  are  used  because,  for  equal  entrance 
area,  the  width  is  such  larger  than  the  height,  and,  therefore, 
any  slight,  hard-to-detect  angle  between  the  probe  width  and 
the  wall  produces  an  upward  shift  in  the  geosetric  center. 
Vibration — Some  experisents  were  sade  to  detersine  if  any 
vibration  was  present  in  the  extresely  thin  probes  during  the 
tests  particularly  when  the  probe  was  near  the  surface.  This 
could  have  an  laportant  effect  on  the  evaluation  of  wall  dis¬ 
tances  especially  for  points  nearest  the  wall.  The  tests  con¬ 
sisted  of  Measuring  the  electric  current  flowing  between  the 


NOLTR  64-99 


probe  and  the  model  wall  when  a  1.5  volt  battery  with  a  10^  ohm 
resistance  in  series  was  connected  between  the  electrically 
insulated  probe  and  the  plate.  The  output  was  recorded  on  a 
sensitive  oscilloscope.  If  vibrations  were  present,  it  was 
expected  that  as  the  probe  approached  the  wall  the  signal 
obtained  would  Jump  intermittently  from  zero.  The  observed 
result  was  a  progressively  decreasing  contact  resistance  as  wall 
contact  was  being  established.  No  significant  vibration  was 
discerned. 

Additional  tests  with  an  unsupported  wire  of  .127  mm  diam¬ 
eter,  about  20  mm  in  length  showed  that  wall  contact  could  be 
repeated  within  +  .025  mm.  This  indicates  the  forces  tending  to 
bend  the  probe  in  the  vicinity  of  the  wall  were  extremely  small. 
Wall  Interference— In  addition  to  the  considerations  regarding 
the  physical  location  of  the  instrument  and  the  accuracy  of 
pressure  transmission  and  recording,  there  are  also  certain 
aerodynamic  effects  that  should  be  considered.  Probably  the 
most  important  of  these  is  that  of  wall-probe  interference. 

That  is,  the  disturbance  to  the  flow,  caused  by  the  probe, 
propagates  into  the  subsonic  flow  near  the  wall  and  there  the 
disturbance  propagates  for  some  distance  upstream.  Under  some 
conditions,  this  upstream  disturbance  may  in  turn  effect  the 
flow  in  front  of  the  probe.  In  such  a  case,  the  presence  of 
the  probe  changes  the  flow  it  is  intended  to  detect^^. 
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Preliainary  lurveys  with  various  size  Pitot  probes  (.89, 
.559,  and  .406  as  in  diaaeter)  were  conducted  with  the  surveys 
extending  through  the  boundary  layer  directly  over  a  static 
pressure  orifice  in  the  plate.  As  night  be  expected,  all  the 
probes  caused  the  static  pressure  to  rise  as  they  approached 
the  wall,  and  the  largest  probe  produced  a  large  disturbance 
which  affected  the  static  pressure  at  greater  distances  froa 
the  wall.  The  aaxiaua  increase  in  wall  pressure  was  also 
larger  than  observed  with  the  two  snaller  probes.  The  latter 
produced  about  the  sane  naxlnun  increase  in  surface  pressure 
(Ap|i/p«  =  25  percent  when  the  probe  center  was  one  dianeter  from 
the  wall).  Since  the  saallest  probe  had  a  prohibitively  long 
response  tine  for  the  low  values  of  Pitot  pressures  occurring 
near  the  wall,  the  .559  an  probe  was  chosen  for  the  nain  test. 

To  deteralne  whether  or  not  the  disturbance  propagating 
on  the  surface  influenced  the  Pitot  aeasureaents,  the  following 
experlaent  was  perforaed.  A  dunay  probe,  l.e.,  a  length  of 
tubing  soldered  closed,  was  located  on  the  aodel  surface  in 
order  to  produce  a  disturbance.  The  Pitot  probe  regularly 
used  for  the  boundary  layer  surveys  nade  a  survey  directly  over 
the  duaay.  Within  the  neasurlng  accuracy,  the  profiles  were 
the  sane  with  or  without  the  duaay  except  when  the  two  probes 
were  in  direct  contact.  It  was  concluded  that,  except  for  a 
region  about  ouo  probe  diaaeter  frow  the  wall,  the  aeasureaents 
are  unaffected  by  wall  interference.  This  justifies  the 
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procedure  of  discounting  sone  of  the  ■eeeurements  nearest  the 
wall  as  has  been  done  by  several  other  investigators^®' 

Turbulent  Fluctuations— Another  possible  source  of  interpreta¬ 
tion  error  involves  the  Pitot  pressure  fluctuations  caused  by 
the  turbulent  nature  of  the  flow.  The  effect  of  the  velocity 
fluctuations  in  increasing  the  Measured  tl«e  average  Pitot  pres¬ 
sure  has  been  considered  for  incompressible  flow  by  «any 
authors®®"®^,  with  respect  to  supersonic  flow,  authors  have 
either  ignored  the  possible  effect  or  corrected  their  neasure- 
■ents  eMploylng  the  incoMpressible  velocity  fluctuation  correction. 
Only  Nothwang®®  considered  the  effects  of  density  and  velocity 
fluctuations  and  then  only  to  explain  qualitatively  the  discrep¬ 
ancy  between  various  Measuring  techniques. 

The  iMportance  of  the  density  or  temperature  fluctuations 
can  be  denonstrated  by  noting  that  the  Pitot  pressure  is 
approxiMately  linearly  dependent  on  the  square  of  the  local 
Mach  nuMber.  The  square  of  the  local  Mach  number  can  be  written: 

=  U Vf  R-  [Tt  -  U Vi c^3  (58) 

In  a  hypersonic  turbulent  boundary  layer,  the  translational 
energy  represented  by  u2/2cp  is  a  large  part  of  the  total  energy, 
represented  by  T^ ,  even  quite  near  the  wall,  and  the  fluctua¬ 
tions  in  the  total  temperature  are  in  general  small.  This 
velocity  fluctuation  affects  the  Mach  number  squared  and  thereby 
the  Pitot  pressure  in  two  ways.  First,  the  pressure  is  affected 
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directly  through  the  square  of  the  velocity  as  in  incompressible 
flow  and  second  through  the  temperature  fluctuations.  In  both 
cases  the  tendency  is  to  increase  the  Pitot  pressure. 

An  analytical  estimate  of  the  Importance  of  turbulent 
fluctuations  was  derived  in  reference  64,  This  calculation 
showed  that  the  error  could  be  quite  large.  For  the  experimental 
situation  considered  in  the  present  investigation,  the  error  in 
the  local  velocity  very  near  the  surface  could  be  as  much  as 
8-12  percent.  This  coeiputed  error  is  based  on  the  maximum 
percentage  turbulent  fluctuations  as  measured  in  incomprsssible 
flow.  This  maximum  la  encountered,  in  practice,  only  when  the 
probe  is  in  direct  contact  with  the  wall,  and  since  the  turbu¬ 
lent  fluctuations  decrease  very  rapidly  away  frrx  the  maximum, 
the  above  errors  are  an  upper  limit  affecting  only  the  points 
nearest  the  wall.  Furthermore,  Klstler®^  has  measured  decreas¬ 
ing  turbulent  fluctuations  with  Increasing  Mach  number. 

Averaging  over  Probe  Inlet  Area — An  estimate  has  been  made  of 
the  possible  error  incurred  by  averaging  the  pressure  over  the 
finite  inlet  area  of  the  probe.  The  measured  pressures  are  con¬ 
sidered  as  corresponding  to  the  locus  of  the  geometric  center  of 
the  probe.  Deviation  from  this  assumption  are  Important  only 
where  the  Pitot  pressure  variation  is  sufficiently  nonlinear, 
as  for  Instance  near  the  surface  of  the  plate,  where  the  varia¬ 
tion  of  Pitot  pressure  with  distance  is  approximately  quadratic. 
Estimates  for  the  conditions  and  the  probe  employed  showed  that 
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the  measured  values  are  at  most  1.4  percent  higher  than  the 
pressure  corresponding  to  the  probe's  geometric  center. 

Slip  Flow— The  conditions  of  the  present  tests  are  such  that 
the  probe  is  not  subject  to  significant  slip  flow  effects, 
because  the  density,  even  near  the  wall,  is  not  sufficiently 
low.  The  local  mean  free  path  is  about  l/30th  of  the  Inside 
diameter  of  the  probe  under  the  worst  condition. 

Conclusion— For  the  circular  entrance  probe  used  and  under  the 
conditions  of  the  investigation,  the  major  errors  are  introduced 
into  the  measurements  only  for  the  last  few  points  before  wall 
contact.  These  particular  points  are  normally  disregarded 
because  of  probe-wall  interference.  The  interference  disappears 
within  one  probe  diameter  from  the  wall.  Turbulent  velocity 
fluctuations  also  affect  the  data  near  the  wall  contact  point. 

An  estimate  shows  that  the  latter  causes  the  measured  values  to 
exceed  the  correct  values.  However,  turbulent  fluctuations  become 
rapidly  less  Important  with  increasing  distance  from  the  surface. 
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POROUS  FLAT  PLATE 


FIG.  I  POROUS  INSERT, 
PLENUM  AND  INSERT 
HOLDER, AND  SURFACE 
PLATE,  IN  THAT  ORDER. 


FIG.  2  PLENUM  AND 
INSERT  HOLDER  SHOWING 
COOLING  FLUID  INLET 
HOLES  AND  INJECTED  AIR 
INLET  SLOT. 


FIG.  3  UNDER  SIDE  OF 
POROUS  INSERT, SHOWING 
STATIC  PRESSURE  TAPS 
AND  THERMOCOUPLE 
STATIONS, THE  COOLING 
TUBES  PROTRUDE  FROM 
EACH  END. 
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FIG.  5  POROUS  FLAT  PLATE  AVERAGE  MASS  FLOW  VS  PRESSURE 
DROP.  AVERAGE  TEMPERATURE  =  20°  C. 
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FIG.  7  EQUILIBRIUM  TEMPERATURE  PROBE  (ALL  DIMENSIONS  IN  cm 


FIG.  9  SURFACE  TEMPERATURE  DISTRIBUTION 
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FIG  14  EFFECT  OF  MASS  TRANSFER  ON  COMPRESSIBLE 

TURBULENT  BOUNDARY  LAYER  MACH  NUMBER  PROFILES. 
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EFFECT  OF  MASS  TRANSFER  ON  COMPRESSIBLE  TURBULENT  BOUNDARY  LAYER 
TOTAL  AND  STATIC  TEMPERATURE  PROFILES 
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DISTANCE  FROM  THE  WALL,  (y  mm) 

FIG.  16  EFFECT  OF  HEAT  TRANSFER  ON  BOUNDARY  LAYER 

TEMPERATURE  PROFILES 
Mo<.=6.7,  Rex=  4X10®,  Cq  =0 
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FIG.  17  BOUNDARY  LAYER  TEMPERATURE 
PROFILE  NEAR  THE  WALL 
Mco  =  6.7,  Rex4X|0®,Ca  =  0 


MASS  TRANSFER  ON  COMPRESSIBLE  TURBULENT  BOUNDARY  LAYER 
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FIG.  25  VARIATION  OF  SKIN  FRICTION  WITH  LENGTH  REYNOLDS  NUMBER 
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MOMENTUM  THICKNESS  REYNOLDS  NUMBER.  Ree^'O 


FIG  26  VARIATION  OF  SKIN  FRICTION  COEFFICIENT  WITH 
MOMENTUM  THICKNESS  REYNOLDS  NUMBER 
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FIG,  28  VARIATION  OF  MOMENTUM  THICKNESS  WITH  DISTANCE 
FROM  THE  LEADING  EDGE  FOR  FOUR  RATES  OF  MASS 
TRANSFER. 
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FIG.  30  REYNOLDS  ANALOGY  WITH  AND  WITHOUT  MASS 
TRANSFER 
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FIG.  38  EFFECT  OF  MACH  NUMBER  ON  VELOCITY  PROFILE  CONSTANT 
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FIG.  44  EFFECT  OF  MASS  TRANSFER  ON  THE  PROFILE  CONSTANT 
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